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Abstract 

We consider a linear Schrodinger equation with a small nonlinear perturbation 
in R^. Assume that the linear Hamiltonian has exactly two bound states and its 
eigenvalues satisfy some resonance condition. We prove that if the initial data is near a 
nonlinear ground state, then the solution approaches to certain nonlinear ground state 
as the time tends to infinity. Furthermore, the difference between the wave function 
solving the nonlinear Schrodinger equation and its asymptotic profile can have two 
different types of decay: 1. The resonance dominated solutions decay as 2. The 

radiation dominated solutions decay at least like 
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1 Introduction 



Let Hq be the Hamiltonian Hq = —A + V — eo with V a smooth locahzed potential and 
Co < the ground state energy to —A + V. Consider the nonhnear Schrodinger equation 

^^t^|J = {-A + V)^P + X\^P\^^|J, ^P{t = 0) = i/jq (1.1) 

where A is a small positive or negative parameter. Our goal is to understand its asymptotic 



evolution as t — > oo. The nonlinear bound states to the Schrodinger equation ( |1 . 1|) are 
solutions to the nonlinear equation 

i-A + V)Q + \\Q\^Q = EQ, (1.2) 

For any nonlinear bound state, ipt = Qe"*^* is a solution to the nonlinear Schrodinger 
equation. We may obtain a family of such bound states from minimizing the corresponding 
nonlinear energy functional: 

m = llm'+lv\<p\' + \x\<p\'dx. 

For each > and A sufficiently small, subject to the constraint ||0||^2 = N, there is a 
unique positive minimizer Q ofTi, which solves (|1.2|) for some E = E{N) and has exponential 
decay as a; — > oo. We call this family nonlinear ground states. 

Instead of N, we can use E as the parameter. From now on we will refer to this continuous 
family as {Qe}e- Let 

He = -A + V-E + XQI . (1.3) 

We have HeQe = 0. Since A is small, the spectral properties of He are similar to those of 
Ho. 

Suppose the initial data of the nonlinear Schrodinger equation ipQ is near some Qe- If 



—A + V has only one bound state, it was proved in W^ that the evolution will eventually 
settle down to some ground state Qeoc with E^o close to E. Suppose now that —A + V has 
multiple bound states, say, two bound states: a ground state 0o with eigenvalue Cq and an 
excited state 0i with eigenvalue Ci, i.e., i/o0i = ^oi4>i where Cqi = ei — Cq. The question is 
whether the evolution with initial data ipo near some Q e will eventually settle down to some 
ground state Qeoc with E^o close to E7 Furthermore, can we characterize the asymptotic 
behavior of the evolution? 

The family of ground states is stable in the sense that if 



inf k'o - Qeg 



j0i 



Il2 



is small, it remains so for all t. Let IHI^^a^oc denote local norm. One expects that this 
difference is actually approaching to zero, i.e., 

lim inf ll^t-gse^^lL^, =0 (1.4) 

In this paper, we shall answer this question positively in the case of two bound states. 
Furthermore, we estimate precisely the rate of relaxation for certain class of initial data. 

We start with a simple question concerning the asymptotic dynamics around a fixed 
ground state profile at t = oo. More precisely, we first fix an E and the corresponding 
nonlinear ground state Qe- Our goal is to analyze the detailed asymptotic behavior of those 
solutions converging to this nonlinear ground state as t ^ oo. Although this problem seems 
to be simple, we found that there are two different types of behaviors for ipt — Qe^^^'- one 
is called resonance dominated solutions ; the other radiation dominated solutions. We first 
explain the Resonance condition on Hq and the meanings of these solutions. Recall that 
the ground state and the excited state to Hq are denoted by 0o and 0i respectively. The 
following two conditions are assumed for this paper. 

Assumption Al: Resonance condition. Let eoi = ei — cq be the spectral gap of the 
ground state. We assume that 2eoi > |eo| so that 2eoi is in the continuum spectrum of Hq. 
Furthermore, for some constant 7o > and all real s sufficiently small. 



1 

Hq — Oi — 2 Cqi — s 



^001 , Im 77 Pc""" 000? ) > 70 > . (1.5) 



Assumption A2: For A sufficiently small and E in a. small neighborhood of cq, the bottom 
of continuum spectrum to —A + V + \Q\, 0, is not a generalized eigenvalue, i.e., not a 
resonance. Also, we assume that V satisfies the assumption in Yajima so that the W^'^ 
estimates for k < 2 for the wave operator Wh holds: for a small a > 0, 

IVVix)] <C {xy^'" , for|a|<2. 

Also, the functions (x ■ V)''V, for k = 0,1, 2, 3, are —A bounded with an — A-bound < 1: 

[X ■ V) V0II2 < ctq II-A0II2 + C II0II2 , ao < 1, A; = 0, 1, 2, 3 . 

Fixed an E and its ground state Qe- Let /^("^ be the operator obtained from linearizing 
the Schrodinger equation ( |1.1|) around the trivial evolution Qe"*^*, i.e., 

£WA; = -i [HEk + XQ\k + k)} , 
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where He is defined in ( [L.3D . It is more convenient to work with operators orthogonal to the 
ground state Q. Let 11 be the projection which ehminates the Q-direction: 

Ilh=h-{coQ,h)Q , co=iQ,Q)-' 

and let X denote its image, 

X = U{L^) . 

Define the operator C acting on the space X = UL^ by 

Ch= -i{Hh + XUQ^U{h + h)} . 

The operator can be written in matrix form 

L_= H, L+ = H + 2AnQ2n. 

With respect to £, we define generalized "eigenspaces" : 

^u{C) := {jp : C^ip = —u'^ip^ = + iv : u,v real, L_L+u = u'^u, L^L_v = u^v^ . 

Since the original Hamiltonian Hq has only two bound states with the ground state projected 
out, there is exactly one value for z/, called k, and dimRE^ = 2. From simple perturbation 
theory, we know that k = cqi + 0(A). We can normalized u and v such that {u,L^u) = k 
and V = K~^L^u. Then {u,v) = 1 and the space is just span^ {u,iv}. 

Define the inner product 

((^,0)) = (Re?/',Re0) + (Im?/',lm0) . 

The space of continuum spectrum of C can be characterized by 

Hc(/:) := : V ^ M^*) for all Jy} , 

where C* is the adjoint w.r.t. the inner product just defined. It is also clear from the 
definition that 

X = E^{C) © MC) . 

Notice that this is not an orthogonal decomposition. Explicitly, the eigenspace Ek[C*) is 
simply spauR {v,iu}. Hence the continuum space is characterized by 

Hc(£) = {wi + iw2 : Wi J- v,W2 -L u} . 



C 
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For any initial data ipo near Q, we can solve 70 and io uniquely such that 

It is more convenient to express ipo in term of Q and R = Re defined by 

Re = OeQe 

We shall see in next section that Re is of order A~^. We now rewrite the initial data as 

ipo = [Q + aoR + ho]e'''', ho ± Q . 
(cf. Lemma 9.1) From the decomposition of X, we can write 

ho = aou + Poiv + r]o, 770 G Hc(£) . 

Denote the norm of ipo by n. Since the evolution is nonlinear, we can always re-scale 
A so that n, the norm of tpo, takes any fixed value. It is more convenient to allow n to be 
a parameter between, say, 1 and 10. Define X\ to be the interval so that the norm of the 
ground state Qe is between 1 and 10. Let zo = ao + iPo. Define the notations 

(x) = Vl {t}^ = e'^ + 2Tt, {t}^ ~ max{£-^t} 

where F is a constant to be specified later on. For the moment, we remark that F is of order 
2A^ times the quantity in (1.5). When the subscript e is understood, we shall drop it. 

Theorem 1.1 (1) [Resonance dominated solutions] There exists small parameters 
Aq and Eq such that for any E El\ and |A| < Aq the following hold. Suppose that 

< l^ol = e <eo, 

||'7o||H2niy2>i{M3) < C'kol^^^ • 
Then we can find a small real number ao = ao{E,ho) such that the solution ilj{t) to the 



Schrddinger equation (|1.1|) with the initial data ipo = Q + clqR + ^0 can be decomposed as 

^{t) = [Q + a{t)R + h{t)] e'®W (1.6) 

with 

a(0) = ao, /i(0) = ho, 6(0) = 0, 
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and a{t),h{t) ^ as t ^ oo in the following sense: Let h{t) = + ri{t) with (^{t) the 
component in Ek(£) and r]{t) the component in Hc(£). Then we have 

\a{t)\ <c{ty^ 

C{t) . . . local, smooth, HCWIL^ ~ {ty'^' , 

7]{t) . . . dispersive wave, h(t)||r4 < , |h(t)L2 < Ct'^ , 

loc 

where {t} = e'"^ + t. Also, e(t) /t -E. 

(2) [Radiation dominated solutions] For any Xoc £ H'^ H ly^'-^(M^) small, there 
exist solutions of the form 

tPix, t) = [Q{x) + a{t)R{x)] e*®W + x{x, t) 

such that 

\a{t)\<Ct-\ M-Ml^ <Ct-'/' 

loc 

and X = X1 + X2, where xi = e'^\ao and \\x2i-,t)\\^2 < o{t-^/^). 

We have used the notation f ^ g for Cg < f < C~^g for some constant C . Thus the decay 
of the excited state for solutions constructed in part (1) is given precisely. Consequently, the 
two types of solutions are completely different. Return to the general question concerning 



the asymptotic mass of the ground state profile at t = 00 (1^). Instead of minimizing the 
L? norm, we prefer to determine the ground state profile by the condition 

(V^t - Qi^W e*®^*) , Qijw) = (1.7) 

We shall prove that there is a unique solution to this equation provided that the initial data 
is near some Qm = Qe^^ (Lemma |9.2|) . Furthermore, we can determine the change between 
E{t) and E-.^. 

Theorem 1.2 Let Aq and eq he given as in the first theorem. For any E[^ G Ix with A < Aq 
the following properties hold: 
(1) Suppose the initial data 

i^O = Qin + «in-Rin + ^in 

satisfies that 
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Then there are solutions E{t) and 0(t) to ( |1.7| ) such that 



limE{t) = E^, \E{t)-E^\<C/{t) 

t— >oo 



Furthermore, if we decompose ijjt as in (pTOl ) with Q = Q^o being the profile at time t = oo 
then we have 

\ait)\ < C {t}-\ hit)=at)+r]{t) , 
C{t) . . . local, smooth, \\Cit)\\L^ < {tV^^^ , 

ri{t) . . . dispersive wave, \\ri{t)\\j^4 < C {t}~^^^ , — Ct^^ . 

hoc 

where {t} = {t}^^_^. 
(2) Suppose that 

<e = \zi^\ < eo, 
||^in||_ff2n^y2,i(iR3) < Ce^^'^, A~"'"|ain| < Ce'^ . 
Then ipQ belongs to the class of resonance dominated solutions with respect to the final profile 

We shall see in next section that i?^; is of order A~^. This explains that has to be 



order A in order Theorem to be correct. Notice that once we have proved the assertion 



(1) of Theorem |1.2| , the second assertion follows from Theorem The proof of Theorem 
pTTI is complicated due to the construction of solutions to the Schrodinger equation with the 
boundary condition of Hq at time t = and that of a at the time t = oo. If we are interested 



only in Theorem |1.2| , we can omit this construction and the proof will be much simpler. We 
feel that this construction may be needed in other contents and so we keep it. 

The resonance solutions for nonlinear Klein-Gordon equations were first proved in an 



important paper [|I^ by Soffer and Weinstein (see also 0). They consider real solutions to 
the nonlinear Klein-Gordon equation 

d^u + B^u = Xu\ ■={-A + V + m^), (1.8) 

with A a small nonzero number. Assume that has only one eigenvector (the ground 
state) (p, B^(j) = Q'^(f), with the resonance condition Q < m < 3Q (and some positivity 
assumption similar to that appears in assumption Al). Rewrite real solutions to equation 



LSI) as u = acj) + Tj with 

a{t) =ReA{t)e'^\ ReA'{t)e'^' 



Then A and rj satisfy the equations 



^ = ^e-™(0,A(a0 + r/)=^) (1.9) 



{d^ + B^)r] = P,X{a4> + r]f (1.10) 



Theorem 1.1 in |12[ states that all solutions decay as 

In particular, the ground state is unstable and will decay as a resonance with rate t~^''*. 
We first remark that the proof in ||T^ has only established the upper bound t~^^^. Fur- 



thermore, an universal lower bound of the form t is in fact incorrect. From the previous 
work of ^, it is clear that radiation dominated solutions decaying much faster than t~^l^ 



exist. Similar to Theorem 1.1, we have two cases: 



1. '/^(O) <^ ^(0) : the dominant term on the right side of ( p..9| ) is . 

2. r7(0) ^ ^(0) : the dominant term is A//^ . 

In case 2, another type of solutions arises, namely, those with decay rate 

We shall sketch a construction of such solutions at the end of section 10. 

Notice that all solutions in |12[ decay as a function of t. Therefore, we can view as a 
study of asymptotic dynamics around vacuum. Although most works concerning asymptotic 
dynamics of nonlinear evolution equations have been concentrated on cases with vacuum 
as the unique profile at t = oo, more interesting and relevant cases are asymptotic dynam- 
ics around solitons (such as the Hartree equations @], see also P). The soliton dynamics 
have extra complications involving translational invariance. The current setting of nonlin- 
ear Schrodinger equations with local potentials eliminates the translational invariance and 
constitutes an useful intermediate step. This simplifies greatly the analysis but preserves a 
key difficulty which we now explain. 

Recall that we need to approximate the wave function by nonlinear ground states for 
all t. Since we aim to show that the error between them decay like we have to track 

the nonlinear ground states with accuracy at least like While the nonlinear ground 

states approximating the wave function can in principle be defined, say, via equation ( |1.4| ) 
or (1.6), neither characterizations are useful unless we know the wave functions precisely. 
Furthermore, even assuming we can track the approximate ground states reasonably well, the 



8 



linearized evolution around these approximate ground states will be based on time- dependent, 
non-self adjoint operators Cf. At this point we would like to mention the approach of ||rT| 
based on perturbation around the unitary evolution e**^", where Hq is the original self-adjoint 
Hamiltonian. While we do not know whether this approach can be extended to the current 
setting by combining with ideas of [|l^] (It was announced in ||12[ that its method can be 
extended to (1.1) as well.), such an approach can be difficult to extend to the Hartree or 
other equations with non-vanishing solitons. The main reason is that these dynamics are 
not perturbation of linear dynamics. We believe that perturbation around the profile at 
t = cxD is a more nature setup. In this approach, at least we do not have to worry about the 
time dependence of approximate ground states in the beginning. But the linearized operator 
C = Coo is still non-self adjoint and it does not commute with the multiplication by i. So 
calculations and estimates based on C are rather complicated. Our first idea is to map this 
operator to a self-adjoint operator by a bounded transformation in Sobolev spaces. This 
map simplifies many calculations and in a sense brings the problem back to the self-adjoint 
case at least as far as estimates are concerned. The price to pay is that the new operator 
involves square root of operators like H. From the standard formula for the square root of 
an operator, we can represent it as an integral of resolvents of H. Thus the linear analysis of 
the non-self-adjoint C is reduced to the analysis of resolvents of H where standard methods 
applied, see section 3. 

The next step is to identify and calculate the leading oscillatory terms of the nonlinear 
systems involving the bound states components and the continuum spectrum components. 
The leading order terms however depend on the relative sizes of these components and thus 
we have two different asymptotic behaviors: the resonance dominated solutions and the 
radiation dominated solutions. Finally we represent the continuum spectrum component in 
terms of the bound states components and this leads to a system of ordinary differential- 
integral equations for the bound states components. This system can be put into a normal 
form and the size of the excited state component can be seen to decay as Notice that 

the phase and the size of the excited state component decay differently. It is thus important 
to isolate the contribution of the phase in the system. Finally we estimate the error terms 
using estimates of the linearized operators. 

The estimates obtained from integrating the equation from t = oo can be viewed as 
uniform bounds for all T, provided that the approximate nonlinear ground states for all T 
are known. On the other hand, to pin down the approximate nonlinear ground states, we 
need to have precise estimates on the wave functions. It is thus nature to consider continuity 
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method by assuming that the approximate nonhnear ground states and various estimates on 
the wave function are known up to time T. We then show that these estimates continue to 
hold up to time T + 6T, with 6T small but fixed, provided that all estimates are re-adjusted 
w.r.t. to the new nonlinear ground state at T + 6T. 

From this outline, it seems that the resonance dominated solutions and the radiation 
dominated solutions occur on equal footing. On the other hand, we believe that the radiation 
dominated solutions are in fact of lower dimension in the space of solutions. A proof is still 
lacking. Although we restrict to the small coupling constant case, it can be replaced by 
spectral assumptions on the operators C with some modifications. The details will appear 
in a future publication. This paper is divided into 10 sections: 



Section 


2. 


The set-up of the problem 


Section 


3. 


Properties of the linearized operator 


Section 


4. 


Main oscillation terms 


Section 


5. 


Estimates on dispersive wave 


Section 


6. 


Excited state equation and normal form 


Section 


7. 


Change of the mass of the ground state 


Section 


8. 


Contraction mapping 


Section 


9. 


Dynamical renormalization of mass 


Section 
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. Radiation dominated solution 



It is a great pleasure to thanks M. Weinstein for explaining to us the beautiful idea in 



the work |T2|] and, in particular, to call our attention to the toy model in [0 which contains 



the basic idea of the resonance decaying in a very illuminating way. 

2 The set-up of the problem 
2.1 Ground state family 

We first review the construction of the ground state family {Qe}e mentioned in Section 1. 
The results we reviewed in this subsection follow from simple perturbation theory and we 
shall not give details. Denote the standard L2 inner product by 

(/, 9)^ f I 9 d'x . 
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Let Q = Qe satisfy 

(-A + - eo)Q + = E'Q 



where E' = E — cq. Let Qe = w(l)o + h with h real and orthogonal to Then h satisfies 

{Hq - E')h + A(w0o + hf = E'w^o, 
We can solve w and h so that 

w^ = 0{X-'E'), \\h\\2 = 0{Xw^) (2.1) 

where we have used that the spectral gap of the operator —A + V — cq is of order one. 
Let R = OeQe- Recall 

lJ'^ = -/\ + V - E + 3XQ^. 
Then by differentiating the equation of Q^; w.r.t E, we have 

L^^^^Re = Qe 

Denote by |0) the ground state to L^°^^ of norm 1. We also have 

Qe = w\0) + O{Xw) 

Hence 

R = (Lf [w\0) + 0{Xw)] = 0{X-'w-') |0) + 0{Xw) = 0{X-'w~^) Q + 0{w) . (2.2) 
2.2 The set-up of the problem 

We consider solutions ipit) to the equation (|1.1| ) The picture is that the solution ipit) can 
be decomposed to two parts: one represents a soliton, the other one the radiation. The 
radiation part will disperse to infinity. The soliton part will converge to a soliton Qncc 
while its L^-norm is changing in time. Hence it is natural to consider solutions of the form 

= {Qeh) + /i(t))e^®^*) , (2.3) 
and study its evolution. If we consider a minimization problem 



inf \\ih — QeG 



j0| 



Il2 



Then = (Qe + h) e*® with Imh Qe, Heh -L Re- Hence we almost have h{t) _L QE{t) in 
our problem, with some small correction. 
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Since this setup would introduce a time-dependent linear operator, we try to find a good 
approximation with a fixed Q = Qe- For this fixed Q, let 

H = -A + V - E + XQ\ 

(we have HQ = 0), and let 11 be the projection which eliminates the Q-direction: 

m = h-{coQ,h)Q , co = {Q,Q)-^ 

Let X denote its image, 

X = n(L2) . 

If we assume 

^ = {Q + A;(t))e*[-^*+^Wl . 

Then the equation for k{t) is 

dtk = C^'^'h - iF{k) - ie{Q + k) , (2.4) 

where 

C^°'^k = -i{Hk + XQ\k + k)} , 

and 

F{k) = XQ{2\k\^ + e) + X\k\^k . 

In view of ( p.3|) , it is natural to assume 

k = aR + h, R = OeQe ■ 

Here a{t) is small and compensates the change of norm of the soliton. Since we would 
like h{t) ± Q for all t. Hence we look for solutions of the form 

ij={Q + aR+h)e'®, {& = -Et + 6(1)) . 

We note that ip can be written in this form if ip is sufficiently close to Q. Specifically, we 
let e*® be the phase of Pqip and choose a so that || -Pq '?/'||j;^2 = \\Q + o,PqR\\l^- Then h is 
obtained hj h = ip e^*® — Q — clR G X. 

We also note, by differentiating the equation of Q w.r.t. E, 

C(°')r = -iQ . 
12 



We substitute k = aR + h into ( ^.41) and obtain 

dth = -dR - aiQ + C^°'^h - iF{k) - i9{Q + aR + h) . 

Since we would like h Q, we choose a and 6 so that the right side of the above equation 
is perpendicular to Q. Using h -L Q and Hh ± Q we obtain the equations for a and 6: 

a = {ciQ,lmF{k)) , ci = {Q,R)-\ 

9 = -{l+coc^'a)-' [a+icoQ,\Q^ih + h)) + {coQ,ReF{k))] , 

where cq = {Q,Q)^^ ■ Then the equation of h is 

dth = Ch + IlFall , 

Ch = -i{Hh + Y{\Q'^Yi{h + h)] , 



F^ii = -i9h-iF{k) - {ciQ, Im F) + ia9 



Ryi , Ru '.= TIR 



Here we have used the equation for a. 
Next we compute F{k) = F{h + aR): 

F{k) = AQ(2|/i|2 + h^) + 2XaQR{2h + h) + SXa^QR^ + X{aR + hf{aR + h) . 

With respect to £, we can decompose X as 

X = E« © (X) , 

where the direct sum is not orthogonal, and 

Ek = {au + [3iv : a, /3 G M} , u,v real, L^u = kv, L^v = ku, (m, v) = 1. 

We write 

h{t) = a{t)u + f3{t)iv + T]{t) . 

Their equations are 

a = up + {v, Re Fall) , 
$ = -Ka + {u, Im Fall) , 

dtf] = Cv+ Pc'-'niFaii) . 

(Note (f , n/) = (f, /), and similarly for u.) The linear parts of the equations for a{t) and 
P{t) form a rotation. To single out the rotation, we study 



z{t) = a{t) + if3{t) = pfi-\ fi{t) 
13 



Here /i ^ captures the rotation part and we expect a slower oscillation in p. Note that z{t) 
can be obtained from h{t) by 



z{t) = {v, Re h{t)) + i {u, Im h{t)) 



The function p{t) satisfies 



{v, Re Fall) +i{u,lmFaii) 

{v, lm(F + 9h) - (ciQ, Im F)Ru) + iiu, - Re(F + 9h) - aORn) 

{v, Im F) - {v, Ru) {ciQ, Im F) + i{u, - Re F) + [{v, Im /i) + - Re h)]e - ic^aO 

({;,ImF) — ReF) + [(vjlm/i) — Re/i) — icza\9 

({;, (F - F)/2i) + -(F + F)/2) + [(^;, (/i - h)/2i) + -(/i + /i)/2) - ic^a\e 



—I 



F) + F) + h) + /i) + caa} ^ 



where C3 = (w, it!n), ■C = — (t", -Rn)ci(5 ^ X, and 

ii+ = -\-v) — 0(1), li- = — v) — 0(A), 
^2+ = + = 0(1), {i_ = 1(1^-5) = 0( A). 

Their orders are derive from the simple facts: {v, Rji) — 0(1), C2 = 0(A) and C3 = (li, it!n) = 
0(1). 

Summarizing, we have 

a = {c,Q,lmF{k)) , Ci = {Q,R)~\ 

9 = -{l + coc^'a)-' [a+ {coQ,XQ\h + h)) + {coQ,ReF{k))] , 
9t77 = £77+P/n(F„„) 



(2.5) 



H ^p 



F) + (m_, F) + {(«+, /i) + (m_, /i) + Caa} 9 



where 



Fall = - iF{h + aR) - (ciQ, Im F) + ia^ 



(2.6) 



F{h + ai?) = Ag(2|/i|=^ + /i^) + 2\aQR{2h + h) + SXa'^QR^ + \{aR + hf{aR + h) . (2.7) 
This is the equation we shall solve for the rest of this paper. 
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Convention: for a double-index a. — [olq^ ai), we denote 

^ 2;"° ^"1, |q;| = cto + cti, [a] = — cto + «i ■ 
For example, z^^^^ = 2;^^^. Hence = /xt^'p"^. In what follows \a\ = 2, = 3, I7I — 4. 

3 Properties of the linearized operator 
3.1 Spectral decomposition 

Recall that, for given E, 

He = -A + V-E + XQI 

and HeQe — 0. Recall 11 be the projection which eliminates the Q-direction: 

Uh^h-(h,Q)Q . 

and the operator C acting on the space X — UL^ by 

Ch= -i{Hh + UXQ'^U{h + h)} . 

We also have L_ — H, — H + 2IIXQ'^II. With respect to £, we can define generalized 
"eigenspaces" : 

Ei,{jC) :— {■^ : C'^ip — —v^tl)^ = + iv : u,v real, L_L^u — i/'^u, L^L^v — v^v^ . 
In our case, the only value for u is k. It is natural to define 

He(£) := {V^ : ± E^(£*) for all v} , 
since it is invariant under C. The orthogonality is defined with respect to the inner product 

{{ipA)) = (ReV',Re0) + (ImV',Im0) , 

and we know 

Ei,{C*) — {v + iu : u,v real, u + iv & E,^{C)} . 
It is also clear from the definition that 

x = eE,(/:)© He(£) . 

V 
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Hence, 



Hc(£) = {01 + i(j)2 : 01 -L V, 02 -L u, for all v, u such that L+L_v — u'^v, L^L^u — v^u\ . 

In our situation, since we only have two simple eigenvalues for Hq and the ground state is 
factored out, we can prove by perturbation argument that there are real u and v such that 



Ek = {au + hiv : a,b & 



Lj^U = KV, L_V = KU, 



1 . 



Since L+L^ and L^L+ are not self-adjoint, it is not convenient to use them to characterize 
Hc(>C). We define two operators 



B = n(L_)^/2n, A = ^BL+B . 

Note that both A and B are self-adjoint and positive on X. With these operators, if we 
define w — k'^^'^Bv, we have w is the normalized eigenvector for the operator A and 

u = i'~^^'^Bw, v = i'^^^B^^w, Aw = Kw, {w,w) = l. 

Hence the continuum spectrum of Hc(>C) can be characterized by 

Hc(>C) = {^jJ = ijji + iip2 '■'4^i-,'4^2 real , _L B~^w^, ip2 -L Bw^ for all v] 
= = V'l + #2 real , S-^i, e Hc(A)} 

= {V' = S01 + iS-V2 : 01, 02 real , 0i, 02 e Hc(A)} . 

The maps B and change the differentiability of the functions. A better way is to use 
A~^/'^B and 74^/^5"^ instead. Let U be the operator 

U(f + ig)^A^^'^B-^f + iA-^/^Bg, for real /,^eX (3.1) 

We will prove in next subsection that these operators are bounded in Sobolev spaces and 
weighted L^-spaces. We have established a one-to-one correspondence between the spectral 
decomposition of X with respect to C and that with respect to A: 

e E,(£) « [/^ e E,(A), 
e Hc(£) < — > Ui/j e Hc(^). 

We have several ways to decompose C as products: 

[B-^ 



1 
-A^ 





A 
->1 



Vo 



u 



J 
u 



B 



.41/25-1 ■ 

A-'/^B 



1 



j_ 

1 



^1 
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We prefer using the last one since it has a simpler form: 

C = U-\-iA)U , 

with U given in 

This decomposition will be especially useful when we integrate integrals of the form 

For this purpose, it is convenient to decompose 

U = U+ + f/_C, U± = A^/'^B-^ ± A-^'^B . (3.2) 

Here C is the conjugation operator, and both U+ and f/_ are self adjoint. Hence [U,i] ^ 
but [f/i,^] = 0. More detailed properties of A and U are collected in next subsection. 

3.2 Lemmas 

Here we collect four lemmas. Recall 

H=-A-E + V + AQ^ A={H^ + IlH^^^XQ^H^/^Iiy/\ 



Also, X is the subspace of all functions in L^(M^) that are orthogonal to Q, and H is the 
orthogonal projection from L^(M^) onto X. 



(3.3) 



Lemma 3.1 (decay estimates for e*^*) For q = 4,8, 
For smooth local functions and sufficiently large (3, we have 



-13 ^-itA 



A-0i-2K 



P/H(a;) 



-/3 



<C{t) 



-6/5 



(3.4) 



L2 



The estimate ( |3.4| ) for A = y/H was proved by Soffer-Weinstein [p!2 
Lemma 3.2 



Qui_ , Im 



A-0i-2K 



^l, Im 



1 



Ho-Oi- 2k 



(Pi) +0(A) >0 . 



(3.5) 



This Lemma is a perturbation result. Notice that if A = then the statement of this lemma 
follows the assumption A.l. Since A is small, by continuity it holds for small A. We define 

1 



T = 2y { Qui , Im 



A-0i-2K 



Pc'^UQul ) > . 



(3.6) 
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Lemma 3.3 (operator U) (a) The operators U and U ^ are hounded operators in W^'T\X 
for k < 2, 1 < p < oo, and in H^'^ fl X for r < 3. (H^'^ is the weighted space with 

ii/Lo„. = ii(x)7iL..; 

(b) The commutator [U, i] is a local operator in the sense 

||[f/,*]n0||^s/.< 0(A) 11011^4 . (3.7) 

We denote the wave operators for C (resp. A and H) by Wc, (resp. Wa and Wh)- 

Lemma 3.4 (wave operators) The wave operators Wc and Wa exist and satisfy W^'^ 
estimates for k <2, 1 < p < oo: (Similar estimates hold for their adjoints.) 



The statement on Wc was proved in following the proof of |T5[. Hence we only need to 
prove the statement on Wa- 

Proof of these lemmas 

We now proceed to prove these lemmas. To simplify the presentation, we will assume 
A > 0. The proof for the case A < is exactly the same. Recall 

H = -A- E + V + XQ\ A={H'' + UH^/^XQ^H^/^uy^''. 

We also denote 

H, = -A-E . 

Note that, if Wa exists, we have the intertwining property that f{A)Pc{A) = WAf{H^)WX 
for suitable functions /. We also have similar property for C. 

Recall X is the space of all functions in L^(R^) that are orthogonal to Q, and 11 is 
the orthogonal projection from L^(M^) onto X. In what follows we will only consider the 
restrictions of H and A on X. Hence we often omit the projection H in the definition of 
A. (It should be noted, however, H^: acts on L^(]R'^).) We denote by if"''" the weighted 
space with norm ||/||j:^o,r = || (3;)'"/||^2. In the remaining of the subsection, when we write 
L^, W^'^, or H^'^ , we often mean their intersection with X: f] X, W^^^ n X, or H^''' n X. 

Assumption on V: We recall our Assumption A2. We assume that is not an eigenvalue 
nor a resonance for —A + V . We also assume that V satisfies the assumption in Yajima ||r3| 
so that the W^'^ estimates for A; < 2 for the wave operator Wh holds: for a small a > 0, 

|V"\/(a;)| < C {xy^'" , for |a| < 2 . 
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Also, the functions {x ■ V)''V, for k — 0,1, 2, 3, are —A bounded with an — A-bound < 1: 

• V)''V<P\\^ < ao \\-A<p\\, + C \\<P\\, , ao < I, k = 0,1,2,3. (3.8) 

By the assumption, the following operators 

H;'/\x-W)'VH-'/', {x-VfVH:\ H:\x-VfV (3.9) 

for A; = 0, 1, 2, 3, are bounded in L^. 

Since Q is the ground state of H with V satisfying the previous assumptions, Q is a 
smooth function with exponential decay at infinity. Hence the above statements on V hold 
also for Q and Q^. Since V + XQ"^ and V have same properties, in what follows we will 
replace V + AQ^ in if by y and write H — -\-V to make the presentation simpler. So it 
should be kept in mind that the potential V in this subsection is in fact V + \Q^. 

For two operators 5" and T, S is said to be T-bounded if ST~^ is a bounded operator. 
If both S and T are self-adjoint, this implies T~^S is also bounded. A deeper result says 
S^/^ is T^/^-bounded, see [RS2]. We say 5" and T are mutually bounded if both ST~^ and 
TS'^ are bounded operators. This is the case if \\{S — T)T~^\\^j^2 ^^2) ~ < 1 for some 9. (It 
implies immediately \\ST-^\\ < 2. Since ||T0|| < + ||(T-5)0|| < ||^0|| +^||T0||, we 

have \\T(j)\\ < C \\S(t)\\, which implies T is ^-bounded.) 

Lemma 3.5 For each A; = |, 1, |, 2, 3, the operators H^, and are mutually bounded. 

Proof. That and are mutually bounded follows from our assumption on V by 
standard argument. To show and are mutually bounded, it suffices to prove the cases 
k — 2 and A; = 3 by the previous remark. We first show \\{A^ — H'^)H''^\\ < 1, which implies 
the case k — 2. 

\\{A^ - H'')H-^\\ = \\H^/^\Q^H^/^H-^\\ < \\H^'''\Q''H-^\\ < \\hI/^\Q''H;^\\ < 1/2. 

The last inequality can be obtained by writing hI'^Q^H;' = H^^'^Q^ + hI'\Q^,H-^] = 
ii-'/'g2 ^ i/-V2[g2^ H,]H-\ and noting [Q^, H,] = AQ^ + 2VQ^ ■ V. 

To prove the case A; = 3, it suffices to prove A^ < CH^ and H*^ < CA^. Note 

{fA^f) = ifA'A'A'f) < ifA'H'A'f) . 

Since A^^H^ + H^I'^XQ'^H^I'^ , we have 

{fA^H'^A^f) < C{fH''H^H^f) + C{f{H'/''XQ^H'/'')H''{H'/''XQ''H'/'')f) 
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where the cross terms are estimated by Schwarz inequahty. To show that the last term is 
bounded by C{fH^f), we shall show that if^/^Q^if"^/^ is bounded in X. Rewrite 

^3/2g2^-5/2 ^ /j3/2^-2g2^^1/2 ^ ^3/2 [g2^ H-^]H~'I^ 

Since [Q^, H"^] is of the form J2\a\<3 ^'^(a;) V", the operators on the right side of the equation 
are bounded in X. This shows < CH^. That < CA^ is proved similarly. Q.E.D. 



Recall the standard formula: 



1 ds 



< a < 1 . 



(3.10) 



The operator T in the above formula will be A^ or H . Hence we also need to estimate 
operators of the form j^jj2- Clearly, for s > 0, 



Lemma 3.6 Let s > 0. The operator 



< 1, 



H 

JTiP 



IS 



I \\(w^'P,W^'P) — 



hounded in W^'^ fl X with 



(3.ii: 



H 



Also, for k = ±1,±2,±3, 



1 



H 



<C{s)-' 



'W'''P,W'''P^ 



H 



<c{sy 



-1/2 



(3.12) 



(3.13) 



Proof. We can rewrite 



H 



- + 



Therefore, to prove statements for j^j^, it suffices to prove the corresponding statements 
for . . We first prove ( |3.12| ) for A; = 0. Let ki denote the eigenvalue of the excited state 
of H and Pi denote the projection onto the corresponding eigenspace. We can write 



H 



St 



X 



Ki + ^ys^ 



■Pi + Wh 



1 



E + yfsi 



H 



where p = —iV and Wh is the wave operator of H. Note E < 0. Since Wh and W^ 
are bounded in W'^'^ for sufficiently nice V, (JlSl), it is sufficient to prove that p2_Ej.^i are 
bounded in W'^'^. However, p2_J"-|-^^ are convolution operators with explicit Green functions: 



C 



\x\ 
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Since |e l^'*^ e±^/si)^/^^ ^ g c\x\{s)'^^^ ^ L^-norms of the Green functions are bounded by 
(s)'^^"^. By Young's inequality we have 

1 



p'^ — E ± ^/si 



[LP, LP) 



which proves (|3.12|) for /c = 0. For k>l and for cj) G W^'^, we have 



H 



W'^'P 



H 



LP 



H 



-1/2 



LP 



W^'P 



This proves ( |3.12| ) for /c > 1. 

For ( 3.131) , we prove the second part. The proof for the first part is similar. For A; > 0, 



since 



[(a;)"' ,H + y/si]= 2V*(V {x)") - (A {xY 



we have 



{xY 



< c 



(V* + 1) 



[x) 



k-l 



(x) 



~k+l 



< C (s)-'/' 



by induction in k. We have the same estimate for (x) , jj}^^ and hence ( |3.13|) holds for 
positive k. The proof for the case A; < is similar. Q.E.D. 

Recall is the weighted space with norm ||/||j:^o,r = || (a^)''/|lL2. 

Lemma 3.7 The operators H'^'^A-^'^, A-'^'^H^'^, R-^'^A^/^ and A^'^H^^/^ are hounded 
operators in W'''^ fl X and H^'^' fl X. 

Proof. By (|3.1CI|) we can write 



^1/2^-1/2 ^ ^1/2 



1 ds 



^1/2 



1 



+ 



s + m s + 

H 



n=0 ^ 



H 



ds 



/•CO TT / TT \ ^ TT 



ds 

7T74 



Since < (s) by Lemma we have 



.-1/2 



poo ^ 
^ ' ' Jo n=0 



-1/2 



ds 
7T74 



< 1 + CA . 
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Similarly 



Also, using (p.l3|) , for r < 3 we have 

||_f/'l/2^ 1/2||^^^^^^^^^^^ _j_ 1/2^1/2||^^^^^^^^^^^^^ ^ -\-CX . 

The above proves that H^^'^A'^^'^ and A~^/'^H-^^'^ are bounded in W'^'^ and H^''^. Indeed, we 
have proved 

11(^)3(^1/2^-1/2 _ 1) + \\{xf{A-'/'H'/' - 1) (x)^ll(^,^^,) < C\ . (3.14) 

We now consider R-^/^A^/^ and A^^H-^/^ Since = AM'^/^ = /^^ ^ 



fi' we 



have 

1 



^-1/2^1/2 ^ ^-1/2(^2 ^ ^l/2;^g2^1/2) T 



S + + m/2XQ2jJl/2 ^3/4 



The main term is Ji. The term J2 is similar to H^^'^A^^^'^, and its integrand has a better 
decay in s for large s. Hence 

||Ag'/2|| <CA||/2|| <cx . 

For the main term /i, 

n=0 



Hence 



Jo „=o ^ "^0 ^ 

Here the norms are taken in (ly^'P, ly^-P) and (if"''", iJ^'*"). Hence we have proved Lemma 



3J. Q.E.D. 

In fact, the last part of the above proof also shows 

\\{xf{H~'/'A'/' - 1) (x)^ll(^,^^,^ + \\{xf{A'/'H-'/' - 1) (x)^ll(^,^^,) < C\ . (3.15) 
The above lemma proves part (a) of Lemma p.3| stating that U and are bounded 



in W'^'P and Moreover, (|J§ and (|37[5| ) mean that f/ - 1 and f/^^ - 1 are "local" 

operators. In particular, they imply part (b) of Lemma |3]^ and that, for any G L^, 



(f/=^i - l)e-**-^-0 ^0 in L^, as t ^ 00 . (3.16) 
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We now prove Lemma We only need to prove the statement on Wa_- Notice 



Wa = lim e^*^e-**^* = lim [/e^^t/^^e^^*^* = lim t/e'^e"^*^* + lim [/e*^(f/-^ - l)e-'*^* 
By ( |3.16|) we have 

Wa = hm Ue'^e'''^* = UWc • 

t^oo 



The bounededness of Wa follows from that of U and Wc- This proves Lemma 3.4. 



We now prove Lemma 3.1. Since 



e"^*^ = WAe-''"'WX Pc^ 



c 1 



the estimate ( p.3[ ) follows from the usual (L^, L'^) estimate for e **^* and the boundedness 
of Wa and P/^ in -spaces. To prove ( |3.4| ), either we prove the bounededness of Wa in 
weighted spaces or we use the Mourre estimate. We will follow the second approach 

and the argument in ||T2 . 



Let a = 2k. We consider intervals A = {a — r,a + r). Let gA{t) = go{{t — a) / r) , where go 
is a fixed smooth function with support in (—2, 2) and go{t) = 1 for |t| < 1. We will consider 
with r small enough. Let D = xp + px, {p = — iV), and the commutators 

ad?) iA) = A, ad^+' (A) = [ ad^ (A) , D] . 

We need to prove the following lemma. 

Lemma 3.8 For A small enough, the Mourre estimate 

gAiA)[zA,D]gA{A)>9gA{Af 

holds for some 6 > 0. Also, gA(A) ad'^(A)g/^(A) are bounded operators in L'^ fork = 0, 1,2,3. 

We will use the following lemma. 
Lemma 3.9 The operators 

are hounded in L? , for k,m = 0, 1, 2, 3. 
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Proof. This is standard and we only sketch the proof. If m is even, we can compute the 
commutator [D'',H"^^'^] exphcitly and estimate 

If m is odd, we write 

^ Jo S + H' ' y/^ 



and proceed as in the case m is even, by using (|3.13|) . Here we have used the formula (|3.1CI| ). 



Q.E.D. 

Proof of Lemma |3.8| 

Let G = A — H and we write A = H + G. Since 

[iA, D] = [H, + V + G,iD] = -A + [V + G,tD] = A-V -G+[V + G, iD] 

and gA{A)AgA{A) > 2egA{Af for some 26 > 0, it suffices to show that, for M = -V + 
[V,iD],—G and [G,D], the operators 

gAiA)MgAiA) = {gA{A)Hl) {H;'MH;') {Hlg^iA)) 

are bounded by gi^{AY and the bound goes to zero when the interval A shrinks to zero. Since 
both gA{A)H^ = {gA{A)A^){A-^H^) and H^gA{A) = {H^A-'^){A^gA{A)) are bounded and 
converges to zero weakly when A shrinks to zero, this will be true if one can show that 
H-^MH-"^ is compact. The case M = —V + [V, iD] is standard and follows from our 
assumption, so we only consider H^^GH^^ and H~'^[G, D]H^'^. 

We proceed to find an explicit form of G. By ( |3.10|) with T = A"^, a = 1/2, we write 



A-' 



ds 



+ 



where 



n=0 ^ ' * 
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By Lemma || Jo II (l^.l^) < /o°° (^)"^^^ ' ^ " {s)'^'^ 3-^'^ ds < C\. Hence 



Since H (a;) ^ and (x) ^ if are compact, from (|3.17| ) ^GH^ ^ is compact. We 
can also write 

h;'gdh;' = {h;'gh'/' (x)} ■ {{x}-' r-'/'dh;'} 



The second operator is bounded by Lemma |3.9| . The first is compact since its terms are of 
the form: H^"^ ■ {x)~^ ■ (bounded operator). Similarly H^^DGH^^ is also compact. Hence 
we conclude the Mourre estimate. 

To show that g/s,{A) a.(^j~i{A)g/^{A) are bounded for /c = 0, 1, 2, 3, we rewrite 

g^{A)B.d],{A)g^{A) = {g^{A)A^) {A-^H^) {H^^ ^dUA)H-') {H^A-^) {A^g^{A)) . 

We only need to show that H^^ ad'^{A)H~^ are bounded since the other terms are bounded 
by Lemma |375| . Recall A = H+G. It is standard to prove that H-^ad'}j{H)H-^ are bounded. 
For H-^ad''^{G)H-\ since it is a sum of terms of the form 



H-^D''GD"'H-\ 



k + m < 3 



it suffices to show that these terms are bounded. By the explicit form ( |3.17|) of G and Lemma 
|3^ , they are indeed bounded. For example, 

^-3^2 1^3/2 ^^^-3 ^^y3 ^-1/2| ^1^-3 

= {H-^D^H""^ (x)-'} Jo{{x)-''H-"^D'H-^] , 
a product of three bounded operators. We conclude that gi:^{A) a.d\){A)g^{A) are bounded 



for k = 0,1,2,3. 



Q.E.D. 



With Lemma |3^ , (cf. the remark in [12], p. 27), the minimal velocity estimate in and 
Theorem 2.4 of ||10l implies 



F{D<et/2) e~'^'gA{A) {D) 



-3/2 



(L2,L2) 



<C{t) 



-5/4 



The same argument in [^] then gives the desired decay estimate ( p.4| ) in Lemma ^TT 



Finally we prove Lemma p.2| . Let ipx = P/^HQu^ and ipo = Pc^° 0o0i- Recall iio = 
-A + V — Cq. We have ipx = ^^0 + 0(A). We write ip\ = ipo + bcpi + r], where t] G Hc(-fio) and 



25 



b,r] = 0(A). Rewrite 



oo 



oo 







Imi {tPx, e-''^^'-^'-^^^i:x) dt (3.18) 
Jo 

poo pt 



'0 ^0 

The main term lies in ( |3.18 ). It is 



Im I j (^0 , e-^*(^«-o*-2«) V-o) rft = ( V^o , Im _ _ ^^ ^c 



which is the desired main term in Lemma 3.2. We want to show that the rest of (B.ISI) and 



( p.l9| ) are integrable and of order 0(A). Recall we write 'j/'a = "^o + + V- Fo'^ term rj 
in by decay estimate we have 

I (V'A , e-^*(^"-°-^'^)r^) I < C {t)-'/' miL^nL^ Wvh^nL^ < C {t)-"' A , (3.20) 

hence this term is integrable. Also, since i^o^i = eoi</>i, 

so we can integrate this oscillation term explicitly. (The boundary term at t = oo vanishes 
due to the decay of e~**(~°*\) We conclude that the rest of ( |3.18| ) are integrable and of order 
0(A). 

For ( |3.19| ), it suffices to show its integrability since AQ^ + G gives the order 0(A). Rewrite 
the last ip\ in ( ^.19| ) as 60i + Pc^° ip\. For the part containing we have 

(^A , e-'(*~^)(^-°'-2'^) (AQ2 + G) e-'^(^»-°'-'")60i) 
= (V^A , e-^*(^-°^-2-) e^-(^-^oi) (AQ' + G) . 

Integration in s gives 

{{A - eoi)- Va , e-^*(^-°^-2'^) {XQ^ + 0) h<P,) . 

Since Cqi lies outside the continuous spectrum of A, the last expression is integrable in t 
following the same argument as ( p. 201 ). For the part containing P(/^° ipx, since (AQ^ + G) is 



a "local" operator in the sense that it sends functions to L^, we have 
I (^A , e-^(*"^)(^-°-2'') iXQ' + G) e-(^o-o.-2«) p^o | < (t - .)-^/^ (.)-^/^ II^aII^hl^ 
which can be integrated in s and t. Hence we have proved Lemma T^. 
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4 Main oscillation terms 



We now identify the main oscillation terms in equation ( p.5|) . We shall use the complex 
amplitude of the excited state, z, as the reference. Recall z(t) = e~'^'^^p(t). Eventually we 
will have p(t) ~ and its oscillation is much smaller than k. The change of mass on 

the direction of the nonlinear ground state is given by a. We will also show that a = O(z^) 
and the order of the dispersive wave 77 (t) is also of order 0(2;^). Assuming these orders, the 
second order term in F, F^'^\ is given explicitly by: 

= AQ(2|CP + C') = ^'0(20) + ^#(11) + ^'0(02) 
where C = zu^ + zu- and 

0(20) = ^Q{u\ + 2m+m_) = \Qul_ + 0(A2), 

0(11) = 2AQ« + ul+ = 0(A), (4.1) 

0(02) = XQiul + 2u+u^) = 0(A2). 

We shall write F*^^^ = z°'(f)a, where a is a double indices (ij) with i+ j = 2 and i,j > 0. 
The repeated indices mean summation. We shall use (3 later on to denote double indices 
summing to 3 and 7 for 4. 

4.1 Main oscillation term in a and F 

We start with identifying the main oscillation terms of a{t). We have fixed the boundary 
condition of a at t = 00 and set a (00) = 0. Thus we can rewrite the equation for a{t) into 
the following equivalent integral equation: 

a{t) = I {ciQ,lmF{k))ds . 

As the oscillation term of order z'^ in F comes from F^'^\ we have up to second order in z'^ 

{ciQ,lmF) ~ 

where 

= (ciQ, AQImC') 

Since ( = zu+ + zu^, we have 

Im ( = Imz (m+ — u_ ) , Im (^^ = (Im ) (n^ — ) . 
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Therefore, we have 

^(2) = Ci IIn^^ Ci = {ciQ, \Q{ul - M^)) = 0(A2) , 
We can integrate A^"^^ by parts to have: 



/ ^(2) ds = Ci Im / ^2 = Im /" (is 

^ oo J oo J oo 



1 



Ci Im { ji ^p^ — iJ, ^2pp ds 

-2iK 



As we shall prove later on, the last integral is higher order term. The first term on the right 
hand side can be written explicitly as 

1 



Ci Im -^—n V = 020(2;^ + z'^) , 
—2zK 



where 



a2o = ^ = ^ {ciQ, Q{ul - u'_)} = 0{\^) , (4.2) 

We shall prove later on that the last term 020 (-^^ -\-z^) is the main oscillatory term in a. 
We denote the rest by 6, i.e., 

a = a2o(-z^ + + 6 (4.3) 

As shall prove a,b,d = 0{z^), but b = 0{z^). In other words, b is the part of a that has 
slower oscillation. From the equation of a, we have the following equation for b: 

b = (ciQ, Im(F - F(2))) - 4 Re a2o-ZA*"^P ■ (4.4) 

Assuming that b and r) are of order z"^, we can decompose F into 

p ^ ^(2) ^ ^(3) ^ ^(3) ^ ^(4) 

where F^^^ and F^^^ denote terms of order 2;^ and z^, respectively, and F^^^ denotes higher 
order terms: 

F(2) = Ag(2|C|' + C') = ^'0(20) + ^^0(11) + ^'0(02) 
F(^) = 2Ag[(C + C)^^'^ + Cn^'^] + AlCrC + 2Aa2o(^' + z')QR{2C + C) 

F(^) = 2\bQR{2C + C) (4.5) 
F(^) = 2AQ[(C + C)^^'^ + C^^'^] + AQ [2\ri\^ + ri"] + 2\aQR{2ri + r;) 
+ 3AaW + A0A;rA;-|C|\] , 
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We can also rewrite the equation of into 

e = C2{z + z) + Fe , (4.6) 

where C2 = — (cqQ, XQ'^u) — 0{X) and 

Fe = - — {coCi^C2a{z + z)+\a + {cqQ, XQ\rj + fj)) + (cqQ, ReF(A;))l | 
1 + CqC^ a 1. L J J 

Hence, since z = iJ,~^p, 



9{t) 



I 2c2Re{z) + Feds = 2c2Rez/{-iK) + f -2c2Re{id~^p)/{-iK) + Feds 
Jo Jo 

^—lmz+ [ -—lm{fi-^p)+Feds (4.7) 

K, Jo 1^ 

4.2 Main oscillation term in 77 

We now identify the main oscillation term in 77. We first rewrite the equation of 77 using the 
operator A as 

dtn ^Cri- Pc^ iOr) + Pc^ UF^ , 

F« = -ieC - iF(k) - [(ciQ, Im F) + iae]Ru ■ (4.8) 

Notice that F" and F*'' differs only by the term iOrj. Observe also that —i9( is not killed by 
P/ since [£, i] 7^ 0. Let 

77* = C/77 . 

Since C = U-\-iA)U and C/P/ = P/ U, we have 

9^77* = -iArj" - Pc^ UieU-^rf + Pc^ C/nF« 

= -ar/^ - lOrf - P/ [[/, ij^tZ-^r;^ + Pc^ UHF^ . 

Let 77 = 6*^77* and use U~^rf — 77, we get 

dtfj = -a?7 + e*^ Pc^ t/HF" - e*^ Pc^ [U, ijOrj . (4.9) 

Hence 

fi(t) = e-*^*77o + e-*^(*-") Pc^ |e^^C/HF« - e^^[C/, 71^77} ds . (4.10) 

Since 77 = U~^e~^^'rj and C/ is bounded in Sobolev spaces, for the purpose of estimation we 
can treat 77 and 77 the same. 
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From the the definition of r/ ( 4.1U ), the integrand is e^^f/IIF^ — e*^[f/, i\6rj. We first identify 
the main term in F^: 



F« = -iOC - iF{k) - [{ciQ, Im F) + iae]Rn 
= iz''(f)l + F^K 



(4.11) 



We have aheady decomposed F into orders in z. To identify the main term of F^, it remains 
to decompose —iO( and (ciQ, ImF)i?n- From the equation of 6 we have 



—i9( = —ic2{z + z){zu+ + zu^) — iFg( 



Also 



(ciQ,ImF)i?n = (ciQ, 



^20 - 902 



)Rnlmz^ + 0(z^). 



Recall the decomposition of F in (^4.5| ) and F'-^^ 



Now we have (|4.11|) with 



= -iFeC - i{F - - [{ciQ, Im(F - F^^^)) + ia9]Rn 



and are defined as follows: 



020 = -(t>20 - C2U+ + |(CiQ, 020 " 0O2)-Rn = "020 " C2M+ + 0(A^ 
011 = -011 - C2(m+ + M-) =0(A), 

C2M_ + i(ci(5, 020 - 0O2)^n = C'(A^) . 



(4.12) 



^02 



^02 



Here we have used 02o = 0{X) = 0ii, 0o2 = O(A^), C2 = 0(A) and -Rn = 0(1). Also note 
that they are all real. 

Recall (U). Since Ui = (f/+ + U^C)i = t{U+ - U^C), We have 

V^e'^Uiz''Y{(t)l = P/e^^i(t/+ - U^C)z'^U(f)l = e'^iz''^^ 



where 



$20 

$02 



{c/+n0«o - fZ-n^J^} , $11 = {(f/+ - u^Wli} , (4.i3) 



-f/_n0«o + f/+n0i 



02 



Hence we can rewrite the integrand in (|4.10| ) as 
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The leading term in rj is 

Jo 

;-**^e*^(^-°') (/i["l(eV)(s)«$a) ds 



itA 



i{A -Oi+ [a]K] 



s=0 



^(2) _ ^-itArje^a 



e*V)(0)^, + (/J) 



where 



^(2) _ e'^z"^^, 



1 



A - Oi + 



(4.14) 



(4.15) 



and (//) is the error from integration by parts, 



(//) = - /'"*e-*(*-^)^ L"A(eV)^a \ ds . 



Also note the sign of Oi is so that the last two terms in ( [1.14| ) decay as t ^ oo. 

We have identify the main oscillation term in fj and we denote the remaining term in by 

(3). 



(4.16) 



Notice that f]^"^^ is not in and this is not an decomposition. It is, however, very useful 
for local behavior as it identifies the local oscillation. 
From ( [4.10| ) we obtain the equation for f]^^^: 



^(■^)(t) = e-'^% - e-^^*(e^^z")(0)^„ - / e 

'o 

ft 



^^(*-^)//i["l^(eV)^4 ds 



+ f e-'^^'"'^^^^e'^UIi{-irfff)ds 
Jo 



(4.17) 



^(3) , ~(3) , ~(3) , ~{3) , ~(3) 

= nil' +r]\' + +r]\' + rill' ■ 

We treat flf^ separately because rffj is a non-local term. Notice that fi^"^^ ^ L^, but is still 
"orthogonal" to the eigenvector of A. 

From the decomposition of fj and the relation t] = U^^e^^^f], we have the corresponding 
decomposition for 77: 



(3)1 
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where 



V 

Summarizing, we have decompose a, F and t] into terms in order of z. The main oscilla- 
tory terms of order z"^ in a is a2o{z^ + z^) and 



with 



a = a2o{z + z ) + b 



b= (cig,Im(F-F(2))) -ARea2oZfi-^p 



(4.18) 



(4.19) 



The nonlinear term F is decomposed into orders in ( [4.5|) with the second order term F^'^^ 
explicitly given. We also rewrite the equation of 6 into 

e = C2{z + z) + Fe , (4.20) 

where C2 = {cqQ, XQ'^u) = 0(A) and 

Fe = ^ , \coC^'c2a{z + z) + \a + {cqQ, XQ^r] + r])) + (cqQ, Re F{k))] ] (4.21) 
1 + CqC^ a I. L J J 

The dispersive wave t] is related to rj by the relation r] = U~^e~^^fj with fj satisfying 

fl{t) = e-'^% + e-*^(*-') Pc^ |e^^f/nF« - e'^U, i]9i]^ ds . (4.22) 
Furthermore, it can be decomposed into t)^^-* + r/*-^^ with 

1 



A - Oi + [a\K 

and r]^^^ satisfies the equation ( [4.17 ). For Theorem LI the boundary conditions are 



a(oo) = = 6(oo) 
< kol < £^0 



ho||//2niy2,i(M3) < \zq 



(4.23) 



We now collect a few properties of the operator U . We can expand U± in order of A as 

;7+ = l + 0(A), f/_ = 0(A). 
Notice that [/c2n+ = 0(A2) since f7M+ is almost orthogonal to He (A). Hence we have 
$2o = -P/n02o + O(A2), $11 = 0(A), $o2 = 0(A2) _ 
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We may decompose 



(4.24) 



then 



where 



V20 = Ulri2o + Ulr]o2, r]n = U ^r]n, r7o2 = f/+?7o2 + f/!^?20 



If we expand 



U+ = 1 + 0{X), U-=0{X), Ul = l + 0{X), U1 = 0{X), 



then we get, by (|]T|) and (|4A3|) , 

-1 



V20 = V20 + 0{X 



A-0i-2K 



P/n02o + 0(A2) , r/n = 0(A) , r/02 = OiX' 



5 Estimates of dispersive wave 



We now estimate solutions to the equations ( |2.5|) with the decompositions into main oscilla- 
tory and higher order terms in section 4. We first need to choose a suitable norm. Define 

{t} = 5-2 + 2n, {t} ~ max{e-^t} . 
and, for 7]^% = Ei=2^f ^ = ^7^^^ - Vi\ 



0<t<T 



M{T):= sup {{tY" + {tf"- Ut)\y + {ty^''" v?\{t) 



7'2 



(5.1) 



Recall that 6, the slow oscillation part of a, is defined via ( ^.18| ) and satisfies the equation 
([4.19|). We assume that the function b is in the following space Bt'- 



Bt = {b{t) : \b{t)\ < D {t}-^ ,0<t<T} 



(5.2) 



where D = 2i?22/r is some constant. The precise form of the constant is not important, it 
merely has to be an order one constant bigger than the true behavior of b which we shall 
derive. The class B is nonempty since it contains the constant function 0. 

Our setting is thus given at the end of section 4 except we now assuming the estimate 
on b (and thus on a too) up to time T. From now on, we fix T. 
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Theorem 5.1 Suppose r],h,a,6,z are solutions to the equations ( ^.51 ). Assuming the esti- 
mate ( |5.2| ) on h, we have 



M{t) < 2 



for all t <T . 



Moreover, if we further assume \zq\ = e > and \\r]o\\ < e^^"^ , then \z{t)\ > c{t} 



Our strategy is to show that M(0) < 3/2 and that M{t) < 3/2 if M{t) < 2. By continuity 
of M{t), this would imply that M(t) < 3/2 for all t < T. So for the rest of this and next 
sections, we shall use freely that M{t) < 2 to prove that M{T) < 3/2. The proof of this 
theorem will be completed after Lemma |6.1| . We now start the proof. 

Due to the presence of the non-local term in F, we first need a global norm estimate 
on rj, which we choose to be ||r7(t)||2,4. Our goal is to prove that 



which agrees with that of free evolution. 

We first recall some basic facts concerning the Schrodinger equation which provide some 
useful feeling on the size of various quantities. Since we shall proceed with iteration scheme, 
our a, z and r] do not solve the Schrodinger equation and we will not use these facts. 

The norm of is uniformly bounded if A is sufficiently small. (It can be proved by 
using the conservation of the Hamiltonian and the Gagliado-Nirenberg inequality.) Hence 
is uniformly bounded since h = ipe^^ — Q. Assuming z{t) is bounded, then ({t) and 
ri{t) are both bounded in H^, uniformly in t, since ( = zu^ + zu^ and rj = h — (. 

Return to a global estimate for rj. Since rj = U~^e~'^^fi and U is bounded in Sobolev 
spaces, for the purpose of estimation we can treat rj and fj the same. Recall [f/, i] is a local 
operator satisfying the estimate at Lemma 3.3. 

We will need the following calculus lemma: 

Lemma 5.2 LetQ <d<l< m. {t} = e^^ + 2rt. 



VmL^<C{t}-'^'log{t} 




Also 




If, instead, d > m > 1 
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Proof. Denote the first integral by (I). If t < £ ^, tlien {t} ~ e ^ and 
Jo 

If t>6-^, then {t} ~ 2Tt and 

(/)</ Ct-^{sY'^ds+ C\t - s\"'C {t}-"^ ds 

Jo Jt/2 

< Ct" ^("^-1) + Ct'-'' {t}-"" ~ C {t}-'' + c {t}'-''-"' 

For the second case, denote the second integral by (II). If t < e^"^, then {t} ~ and 

(//) ~ [\t~ s)-'' e^"" ds<e^'^-C {t}"™ . 
Jo 

If t > then {t} ~ T {t) and 

(jj) ~ [\t- s)"^ r-"" (s)-"^ ds < cr-'" (t)-'" ~ c {t}-"" . 

Jo 

We conclude the lemma. Q.E.D. 
5.1 Estimates in and 



Lemma 5.3 Suppose that rj is given by equation ( [4.10|) and recall t] = U ^^fj. Assuming 
the estimate ( ^.2|) on b and M{T) < 2, we have 



Ut)\\L^<c{t}-'/Uog{t} . 

Moreover, we have 

Proof. From the defining equation of fj, we have 

\\m\y < cii^oii^inLV3 {tr'^'+J^c\t-s\-'/'{\\Fi{s)\\^,,, + \e\ imw,] ds . 

From the definitions of and we have 

||i^«(^)IL4/3 + 1^1 \mh < c\e\ \mh/s + \\F{s)\\^,/, + c\ae\ , 

From the Holder inequality and the definition of ||-F(s)||j^4/3 we can bound ||-^(s)||^4/3 by 
11^(^)11x4/3 < C {\\h\\l + \a\ \\h\\,, + \a\'+ \\h^^,,, + \af) . 
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Since \\h^ 



\h\\j^4 and 



ms)\\,,<\\a^ 



Il4 



11^(^)11 



L4 ' 



we have from the assumption M < 2 that 

\\his)\\,, < C{s}-'/' + C{s}-'/Uog{s} <C{s}-'/' , 
Therefore we have ||F(s)||j^4/3 < C {s}^^ and thus 

\\FKs)\\,.;,<c{sy' . 

Since (t) ""^ < e^"^ {t} ^ and ||^7(t)||2,4 ~ 11^(^)11^4, we conclude 

\m\\L^<Ce''^'\\Vo\\mnL^rAt}-'^'+ f \t - s^'/'C {s}-' ds < C {ty'/' log {t} 

Jo 



Here we have used Lemma |5.2| in the last integration. 
We now bound ||^7(t)||^2. Since 

d 



— iv, v) = Re{v, dtv) = Re{v, e'' Pe^ UUF^ - P,^ [f/, i\eri) 



we have 



< C {t}-^/' log {t} C {t}-' = C {t}-'/' log {t] 



In the second inequality we use our previous estimates. Hence 



7/4 ■ 



and we conclude ||?7(t)||^2 < s^^'^-, and so is ||?7 



L2- 



Q.E.D. 



5.2 Local decay of 77^^^ 

Recall r/^^-* = U^^e^^^f]^^^ and fj^^^ satisfies the equation ( [4.17| ). We want to show that r/*-^^ is 
smaller than fj^'^^ locally. Define the local norm by 



r2 



1,2 



for a fixed sufficiently large /?o > 0, which will become clear later on. 
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Lemma 5.4 Assuming the estimate ( |5.2| ) on b and M{T) < 2 for all T , we have 



~(3) 
^2-5 



(3) 
^2-5 



/n particular, for a local function (f) we have 



loc 



\{^,v^'^)\<c{t}-'-^/' , 
^,\v\' + \^\')\<c\\v\\% +CM^. Ml<c{ty 



(5.3) 



(5.4) 
(5.5) 



Proof. We first estimate rjj appearing on the right side of the equation for rj^^^ (|4.17| ) 
Note we did not include rjl in the above Lemma: 



~(3) 



loc 



-9/8 



This term would be bounded by e'^^'^ {t}^^^^ if we assumed \\rio\\ < e^^" . However, since we 
only assume ||?7o|| < this term needs to be treated separately. 

For fj^^ and fi^\ the two terms involving rja-, {rja ^ L^), from the definition and the 



estimates Lemma 3.1 on A, we have 



~(3) 
^2 



loc 



<C{t} 



-9/8 



~(3) 



<C j\t- s)"'/' {s}"'/' ds<C {t} 



-9/8 



To estimate , we recall the definition of F"" and rewrite 



(ciQ,Im(F-F(2))) + za^ 



Since M < 2, we have |^-C2(2 + z)| < \a\ < C {sY\ (x)* (F - F^^) - r/^T^) 

C and |(ciQ,Im(F — F(^)))| < {s}^^^^. Hence, by Lemma we have 



< 



L2 



~(3) 
?7l 



loc U 



< I {t- s)"^/' (x)^" (F«» + zt^'t]) (s) 

< / {t - s)-''/^ C {s}-^^^ ds 
Jo 

<C{t}-^^^ . 



L2 



ds 
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We now consider fj^\ Denote the integrand e"*"^*^*"*-* Pc^ e*^f/n(— zr^^f]) by X{t,s) and 
decompose the time interval into {0,t — a) and (t — a,t). From the triangle inequality, we 
have 



~(3) 
^5 



< /* \\X{t,s)ds\\^,^^+ f \\X{t,s)ds\\ 



^loc JO 



t-a 



"^loc 



We can bound the L\^^ norm by either or norm. Thus 



~(3) 
^5 



"^loc 



<C f ^ \\X{t,s)\\^,ds + C f \\X{t,s)\yds 

Jo Jt-a 



From Lemma 3.1, we have 



and 



11^(^,^)11 



11^(^,^)11 



L4 



\t-s\ 



\n^\\ 



From the Holder inequality and the global estimate of r] in Lemma |5]^ we have 

5/2 



h'(^)L4/3< ii^iii.< ({^r'/'in(2 +.))'< c{.r^ . 



7/4 



\it>2e ^, then {t} ~ t. We choose a = t/2 and thus 



~(3) 
^5 



"^loc 



< C 



t/2 



it-"7F^^'^"''^'+^y,/,it-.r 



^ C{sY^ds 



nt/2 pt -, 

<C C {t}-'" {s}-'" ds + C j^^^ Y^^.C {t}-' ds 



< c {t}-'/' + c {ty^ t'/^ < ce'i^ {t} 



-9/8 



On the other hand, if t < 2e ^, then {t} ~ £ ^ and we choose a = t to get 



~(3) 
^5 



< c 



-Cie-T ds 



lo 



Combining two cases, we conclude 



~(3) 
^5 



"^loc 



< Ce^/^ {t} 



-9/8 



The lemma follows from all the estimates on r/ -^'', j = 1 



Q.E.D. 
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6 Excited state equation and normal form 
6.1 Excited state equation 



Return to the basic equation (|2.5|) and recall the equation for p: 

F) + {u_,F) + {{u+,h) + {u_,h) + C3a}e\ , (6.1) 
where C3 = {u,Rn), Ci = {Q,R)~^. Recall that 

z,p = 0{z), a,r] = 0{z'^), r/(^) = O(z^) . 
We now expand the right hand side of the equation for p into terms in order of z: 

p = fi {caz'^ + df3z'^ + dibz + d2bz + P^^^ } (6.2) 

The coefficients will be computed later on and their properties are summarized in the fol- 
lowing lemma. The proof of this lemma is just straightforward computation and the reader 
can check it rather easily. 



Lemma 6.1 We can rewrite the equation of p into the form ( |6.2| ) such that the coefficients 
di, d2 and all Cq, are purely imaginary. Moreover, R.ed2i = — F + 0{X^), with 



r = 2X' i^Qui, Im Pc ng< J > O . (6.3) 

Proof. There are two parts in p equation: («+, F) + ('U_, F) and {(m+, h) + («_, h) + 030} 9. 
We first consider the second part: 

{(m+, h) + (m_, h) + csfl} 9 

= {(«+,C + r7(')) + (u_,C + #) + C3a}- 

■ {c2(z + z) + (coQ, AQ2(r/(2) + r^(^)) + « + (cqQ, ReF^^))} + 
= {{u+,C) + {u^,0) ■C2{z + z) 
+ [{u+, 7^(2)) + (u_, 7^(2)) + csa] ■C2{z + z) 

+ ((w+, C) + («-, 0) ■ [a + (coQ, Ag2(r/(2) + ff^^))) + (coQ, Re F^^))] 

HerePg^"* and P'^^ denote the remaining error terms. We first observe that line 4, lines 5-6, 

,(4) 



and P2 are of orders 0{z^\ 0{z^^ and 0{z'^\ respectively. Hence line 4 contributes to Caz: 



line 5-6 to dnz^ and dibz + d2bz, and Pi^'' to P^^\ We also observe that the coefficients in 



39 



line 4 arc all real. Since there is a —iji factor in front, their contribution to c„ are purely 
imaginary. Next we observe that the coefficients in lines 5-6 are real except those involving 
?7^^) . The terms with lowest A-order are: 

{u^,ri^'^)-C2{z + z) + K,C)-(cog,AQ2(r;(2) + ry(2))) . 

The first part is of order A-'^ since = {P^u+ + 0{\), Pc"^ r]^^^ + O {X'^)) , with Pc^r/^^) = 

0(A). The second part has order A^ terms, but not on z^z. We conclude that the contribution 
to Rc(i2i from the second part of p-equation is of order A^. Finally wc observe that line 5-6 
give a term di^2tiz + d2,2bz , with real di^2 — 0(1) and o?2,2 = 0{X). Since there is a —i/i factor 
in front, their contribution to di and ^2 are purely imaginary. 

Now we look at the first part of the p equation. The contribution to CaZ"' is from F^^^ : 

-i{{u+,z''(j)a) + {U-,z"(f)a)} ■ 

Clearly all coefficients of 2;" arc purely imaginary. Together with the analysis of second part 
of p-equation, wc know Ca are purely imaginary. 
The contribution to dibz + d2bz is from F^^^: 

-i {iu+, 2XbQR{2C + + {u-, 2XbQR{C + 2C)} 

with C = ZU+ + ZU-. Hence all coefficients of bz and bz are purely imaginary. Together with 
the analysis of second part of p equation, we know di and d2 are purely imaginary. 
The contribution to P'-^-* is from F^^^: 

p(^) = -^ [(..+, f(^)) + (s_,fW)] , 

and we have P(^) = ^ - iP^'^l 

The contribution to djsz^ is from F^^^: 

-i[(u^,F('^) + (u_,F^'^. 

We only consider ^21- A coefficient in F^^^ has to have imaginary part in order to have a real 
part contribution to d2i. Hence only the first group of terms in F^^\ 2XQ[{( + C)^^^\ has 
contribution to Red2i- Also, when we decompose 77^^^ = z'^rja, we can disregard 77(11) since 
it is real. Now, 77(20) = 0(A), 77(02) = O(A^), u+ — 0{1) and — 0(A), hence the main part 
of {Red2i)z'^z is in 

-i{u+, 2XQ{zu+)z^r](2Q)) 
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Summarizing our efforts, we have 



Rerfai = Im(M+,2AQu+r7(2o)) + 0(A3) = lm{2XQul,r](^20)) + 0{X^) 



{2\Qul,lm 



-1 



i2\Qul,lm 



A-0i-2K 



P/n02o) + O(A3) 

p/nAgO + o(A^ 



Q.E.D. 



6.2 Normal form 



From Lemma 13.21, we have F > 0. 



Lemma 6.2 We can rewrite the equation (|6.2| ) of p into a normal form; 

4 = S2i\q\'^q + dibq + g . (6.4) 
where q is a perturbation of p given in the proof. The coefficient 621 satisfying the relation 



Re ^21 = Re d2i 



(6.5) 



If we assume the estimate ([57 
satisfies the bound 



on b and M{T) < 2, then the error term g(t) given by ( |6.13| 



m\<c,{t}-'/'-'' 



(6.6) 



for some constant Ci. Furthermore, there is a positive constant a such that \q{t)\ is bounded 
by 



1 - a) {ty'/' < \q{t)\ < {1 + a) {t}' 



-1/2 



and hence 



1 - 2a) {t}-'/^ < \zit)\ < (1 + 2a) {t} 



-1/2 



(6.7) 



Proof of Theorem 5.1 From Lemma 6.1, we have 



{tV^^\z{t)\ < {l + 2a) 



From Lemma 15.3 and Lemma 5.4 we can bound 



{t} 



3/4-0- 
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loc 



by Ce'^l'^. Since e is small, we have proved that M(T) < 3/2 and this concludes Theorem 
§3. Q.E.D. 



We now prove Lemma |0 
Proof. We have 

i> = [co^z"" + dpz'^ + dibz + d2bz + P^^^] 
and we want to obtain the normal form ( |6.4D . We will repeatedly use the following formula: 

at \ IKITL J tKTTl 

where, if a = (aocn), {\a\ = ao + ai = 2, 3, 4 ■ ■ ■ ) 

U{z) = {ao + aiC){p-'p) (6.9) 



+ aiC)z~^ [caZ° + dpz^ + dihz + d2hz + P^^^] 



and C denotes the conjugation operator. The formula is equivalent to integration by parts. 



We first remove Cq,^;". Let 



+ [a\) 



Since [a] is even, 1 + [a] 7^ 0. By ( |678| 



P1=P- CaHZ - — — /iZ fa{z) 



Decomposing fa{z) into two parts, we can write 



dtz'^ = z — ° , iiz^iao + aiC)^; ^Caz" 

1^ m(l + [a]) ^ ' 

r"(ao + ttiC)z"^ [rf^z^ + dihz + ^2^^ + P^^^] 



+ \a\) 



and we get 

Pi = SfSlJ-z^ + dinbz + d2nbz + nP^'^^ + gi 

with Sf3 = d/3 + dp. Since c?^ are purely imaginary, due to that Cq, are purely imaginary, we 
have the key relation 

Re5f3 = Redf3 (6.10) 

Next we remove d2f^bz. Let 

fid2bz 
P2=Pi- -TT- — • 
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We have 

V2=V\- [id2hz - -^{bp + bp) 

= fi6pz^ + difibz + {fiP^^^ + gi + g2) , 

where 

92 = ~^^(^^ + ^^) • 
Now we deal with Sj^z^^ terms. Let 

E6i3fIZ^ 

Note 1 + [/?] ^ for /5 ^ (21). We have 

P3=P2- fJ'Spz^ + 93 = S2inz'^z + ndibz + (/iP^^^ +91+92 + 93) , 

with 

^ + \f3]) dt^^ ^' 

Finally, since rj^^ = f/^^e~*^ e~*'^*f/?7o is larger than 772-5 when time t is of order 1, we 
need to extract terms of order 0{zr]^^) from fiP'^^K Recall 

^(3) ^ g-iAt~ ^ e-'^'Ur]o, = U~^e~''^f = U~^e-'' e-'^'Ur]o . (6.11) 



Also recall (|6.1|) and ( |6.2| ). In F we have a term 2XQ{{( + C)^^'* + CVi^'^)^ in h a term rj^^ 
itself. Hence in jJiP^^\ terms of order 0{zriQ) are exactly 

P,,„ ^ - ^^(u+, 2AQ((C + CW? + C#^)) - ^y^C(n_, 2AQ((C + CW? + C#)) 
- z/i |(u+,?7f^) + C(u_,r7f^)|c2(z + ^) . 

Clearly these terms can be summed in the form 

/i(# + rif^) + Cix-\z(j) + Z(j), r/f ^) (6.12) 

Here (f) stand for different local smooth functions. Moreover, if we write U-^ = U3 + U4C 
with U3, Ui = 1{BA-^/^ ± p-Mi/2), both self-adjoint, by (|6Tl| ) and 

U-\zf + zg) = z{Usf + U,g) + z{U,g + Uj) 
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(cf. ( |7.2| )), the above is equal to 
The first term can be written as 
with 



We can write the second term similarly as 

at 

Thus we have 

Pzm =PA + 94, Pi ■= Pa,1 + Cp4,2, 94: ■= (fi'4,1 + C5f4,2) ■ 

Observe that ^_q'^_2^ is not in L^, hence we expect slower decay for p^ and g4. Specifically, 
|P4| < ||r?o|| {t}-'^' {t)-'-'' , \q,\ < WmW {t}-' {t)-'-" . 

Now we let 

q = P3-P4 ■ 

We have 

q = 52i^iz^z + ^dihz + (/iP^^^ + 9i+ 92 + 93) - Pz^o + 94 
= S2i\p\^p + dibp + ((/iP(^) - P,^J + ^1 + ^2 + ^3 + 94) 
= S2i\q\'^q + dibq + g , 

where 

g = ((^P(^) - P,^J +gi + g2 + g3 + 94) + 52i(bpp - \q\\) + - g) • (6.13) 

Hence we have arrived at the normal form. (Note that, although we kept p and h in the 
definition of (72, it should be replaced by the corresponding equations (|6.2| ) and (|4.19| ).) Also 
note we have 

^ fid2bz Spjjiz^ 



q = P ^ ~T^ — r~n A^-^ 



Finally, we can check the size of g satisfying the bound ( |6.6D . To conclude this lemma, 
it remains to prove the estimate on q. But this follows from the next lemma. Q.E.D. 
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6.3 Decay and continuity estimates 

In this subsection we present some calculus lemmas which deal with the decay of q{t) and 
its continuity on the error term g{t). We will write 

q(t) = p(t) e-W , 

where p — \q\ and cu is the phase of q. Recall 

{t} = e~^ + 2Tt, {t} ^ max{e~^,t} . 

Before we proceed with proof, we first explain some simple facts of an ordinary differential 
equation. 

Example. We consider real functions r{t) > which solve 

We have the following facts. 

a. All solutions r{t) satisfy 

r{t) <{C + 2n)-^/' with r(0) = C"^/' . 

b. There is a number ri such that if r(0) > ri, then r{t) ~ (C + 2ri)~^/^. 

c. There is a unique global solution ro{t) such that 

ro{t) ~ t~'^ as t ^ oo. 

d. If r(0) < ro(0), then r{t) = in finite time. 

e. If r(0) > ro(0), then 

poo 

/ r(s)^ds = oo. 
Jo 

Lemma 6.3 Let Eq > be small. Suppose p{t) satisfy 

where \g{t)\ < eq {t)'^^^'" , (x>0. We can find < pi < p2 , ( depending on Eq ), such that 

(a) Ifp, < p(0) < p2, then \p{t)\ ~ C {ty^'\ 

(b) IfO< p(0) < p2, then \p{t)\ <C{ty^'^. 



45 



The above example shows that the conclusion of (a) cannot be expected to hold if p(0) 
is too small. This is the main reason that there are two type of asymptotic behavior: the 
resonance dominated solutions given by case (a) and the radiation dominated ones by (b). 
The proof of this lemma is elementary and similar to the next lemma and thus we omit it. 

Lemma 6.4 Let F > 0, cr > 0, and Ci > be given constants, independent of Sq. Suppose 
a positive function p{t) satisfies 

p^-rp^ + g{t), , (6.15) 

(a) Suppose p(0) = e with < e < Eq and 

\g{t)\<Ci{t}-'/'-% {t} = e-' + 2rt. 

Then there is a constant m = m{eo) > 1 such that 

m-^ {t}-^/^ < pit) < m{t}-^^^ . 

Moreover, m{eo) — > 1"^ as Eq — > 0+. 

(b) Suppose that p(0) < Eq and 

\m\<c,{t}^'/'-% {Oo^^o' + 2n. 

Then we have for some constant C 

Pit) < c {t}-'/' . 

Proof. We first prove part (a). Let p+ = m{t}~^^^ and p_ = m'^ with m > 1 to 

be determined. We have p+(0) > p(0) > p_(0). Moreover, 

p+ = -Lm- V+ = -^pI + r(l - m-^)pl > -Tpl + g , 

if 

r(l - m-')m' {t}-'/' > Ci {t}-'/'-" . 

Also 

p- = -Tm^pt = -Vpl_ -T{rn^ - l)pl_ < -Tpl + g , 

if 

r{m' - l)m-' {t}-'/' > Ci {0~'/'~" . 
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Since {t} " < e^'^ , both inequalities hold if 

r(l - m"^)m^ V{m^ - l)m-^ > Ci el"" . 

This is true for m > 1 arbitrarily close to 1 by choosing Sq sufficiently small. By comparison, 
we have p-{t) < p{t) < p+{t) for all t. 

To prove part (b), note we can still define p_|_ as a comparison function. Since p+(0) > 
p(0), the above argument still holds. Q.E.D. 

The following lemma estimates the continuity of p in the error term 'g{t). Here we consider 
the case p(0) = £ > only. 

Lemma 6.5 Let F > 0, o" > 0, and Ci > &e given constants, independent of Sq. Suppose 
two positive function pi{t) and P2{t) satisfy 

Pi^-rpf + Qiit) (i = i,2), pi(o)=p2(o)=£, 

with < s < Eq and 

mt)\ < C, {t}-"^-'' , \5g{t)\ < So {t}-'^'-" , 
where Sg — g2 — gi and Sq > is a small number. Then we have 

\p2{t) - pi{t)\ <Soe^{ra)-' {t}-^'+^^/' . 
Note that the decay rate is improved. 
Proof. Let r = 6p = p2 — pi- Then r satisfies 

r = -Gr + Sg, r(0) =0, G = r(3p^ + 3pr + r^) . 

We have 

r{t)^ e-^^^^^'^'^^6g{s)ds. (6.16) 
Jo 

Note \r{t)\ < p+{t) - p-{t), hence 

G(t) > r[3pl - 3p+{p+ - p-)] = 3r[m-^ - m(m - m"^)] {t}'^ > 2r {t}'^ 
if m is sufficiently close to 1. Hence 

G{r)dr > 2r {r}"^ dr = 2r ^\2r)-^ [In {r}]U > ^ (In {t} - In {.}) 
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and 

g-/;G(r)dr < |^|-(l+a)/2 _ |^|{(l+<x)/2 _ 

By ( |6.16|) we have 

\r{t)\ < f {t}-^'^''^/'{s}^'+^^/'S,{s}-'/'~'' ds 
Jo 

Jo 

Q.E.D. 

7 Change of mass 

For given Qe and h^, we want to find an a(t) = a{E, Hq; t) which satisfies 
a(cx)) = 0, d= {ciQ,lmF{k)), \a{t)\ < C X {t}'^ . 

7.1 The main oscillatory part of a{t) 

We use the following equivalent integral equation for a(t): 

a{t) = I {ciQ,lmF{k))ds . 

J oo 

Note k = aR + h = aR + C + ^- We want to perform several integrations by parts so that 
we get the form 

a{t) = 0{t-^)+ I 0{t-^) ds = 0{t-^) . 

J oo 

As in Section 4, we decompose t] and a as 

= 7^(2) +r^(3)^ a = a2Q{z^ + z^) + h. 
Recall k = aR + C + ^7 ^''^^ 

lmF{k) = lm{XQh'^ + 2XQRah + X{aR + h){aR + h)h} . 

We also denote 

/i-ip = Caz"" + dpz^ + dihz + d2bz + P^^^ = c^^'" + p(=^'^) 
Using the decomposition of F{k), we have 

(ciQ, Im F) = + A^'^^ + A^^") + A^''^ + A^^^ + A^^^ , 
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where 



= (ciQ, AQImC') 
^(^«) = (ciQ, lm2XQRaC) 
A^'^^^iciQ, Im2XQCr]) 
A^^'^ = {c,Q, ImAKPC) 
^(4) + ^(5) = (ciQ, Im {XQt]^ + 2XQRar] + A [\k\^k - \C\\] }) , 

A(^) = (^CiQ, Im { AQ(77(2))2 + 2XQRari^^^ + X [2\CWaR + rj^^^) + C{aR + ^)] }) 
A(^) = (ciQ, Im{Ag[2r^(2)r^(^) + {r,^'^f]+2XQRar,^'^}) 

+ (ciQ, Im {a [2|Crr;(3) + C'^^ + fC + 2|£|2C + }) , 

Since ( — zu+ + zu^, hence 



aR + T) . 



lm( — Imz {u^ — ImC — {Imz ) {u, — u_ 



Therefore 



^(2) = Ci Imz^ Ci = (ciQ, Ag(M^ - u^)) = 0{X^) , 
A^""^ = C2 a Im ^, C2 = (ciQ, 2Agi?(M+ - «_)) = 0(A) . 

First we integrate A^'^^ : 

j A^^) ds = Ci Im /" = Ci Im / //"^p^ ds 

= Ci Im -i- j/i-y - J l^-^2ppds^ 

= ^2Re|^2-2 J zfi-^pds^ 

= a2o{z'^ + z'^) - 4a2o Re f z [caz"' + d^z^ + d^hz + ciaft^ + P''^^] ds 

= 020 (^^ + + /" ^2,3 + ^2,4 + ^2,5 , 
J 00 



where 



a2o = ^ = ^ (ciQ, Q{ul - ul)) = 0(A2) , 
^2,3 = -4a2oRe2;c„2;° , 
^42,4 = — 4a2o Re \df}zz^ + d^bz'^] 
A2,5 = -4a2oRe^p(^) , 
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Note here we have used lie zd2bz = 0. 

As mentioned in Section 4, the term a2o(-2^ + is the main oscillatory part of a. We 
denote the remaining part as b and hence 

a = a2Q{z^ + z"^) + h . 

Although h is of the same order its oscillation is slower since the right side of its equation 

h= (cig,Im(F-F(2))) -4a2oRe;z/i^^j9 



which we obtained from the above computation, is of order O(z^), (cf. ( |4.19| )). 



Another reason for using h{t) is that it is better suited than a{t) for the iteration scheme 
in Section 8. We have 

8a{t) ~ 5o {t}-' log {t} , 5h{t) ~ 5o {t}"' . 
This difference arizes because it is harder to estimate the continuity of the phase. 

7.2 Integration of 0{z^) terms 
we next integrate 

J oo 

= C2lm-L|az-£;.-l(ap)} 

= ^ 2 Re ja^ - y z{ciQ, ImF) + a(c„z° + df^z'^ + dihz + d2hz + p(^))| 

J oo 

c,a = C2 /{2k) = {ciQ, 2XQR{u+ - u.)) /{2k) = 0(A) 
A,a,3 = -c,a{z + z) (ciQJmF^^)) 

A,a,4 = -Czaiz + z) {ciQ , lm{2 XQ RaC + 2XQCv^^^ + A|CI'C)) - 2c,,aRec„2" 
A,a,5 = -Cza{z + z) {ciQ , ImF^^^) - 2c,aaRe{dpz'' + d^bz + d2bz + P(^)) 

We now integrate A'^'^'K Recall U = U+ + U^C, see We will also write f/"^ = 

f/3 + U^C, with U-i,U^ = \{BA-^l^ ± B-^A^I^), both self-adjoint. We will also use the 
following formulas: 

Re j dxf{Cg) = j dx{Cf)g, Im j dxf{Cg) = - Im j dx{Cf)g, (7.1) 
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where 



and (using U = U+ + U_C) 

U{zf + zg) = ziU+f + U^-g) + z{U+g + U.f) 
Recall 7] = U-^T]" = U-^e-'^fj. Thus, by Q, 

= {ciQ,lm2XQCri) = j dx2ci\Q\ri 

= Im / dx2ci\Q'^{zu+ + zu^) {U3 + 1/40)7]" 



where 



(7.2) 



Hence, by (|7j), (think 



Im J dx [(f/3 - t/4C)(#i + #2)] 

01 = 2ciA(5^M+, 01 = 2ciA(5^M-. 
= f/3 and f/_ = -t/4) 
A^"") =lm [ dx (#3 + #4)^/'' , 



with 



^3 = U-sCpi - f/402, 04 = t^302 - ^401- 



We rewrite 



fir) 



where / = e'^^'^fj = "/^o + Jq ^^'^^ ■ ■ ■ ds. The reason we work with /, instead of fj, is that those 
terms in drf of same order {z^) are explicit and do not involve differentiation, compared to 
those in drTj. Now we have 



= Im / (03, zr/^) + (04, ZT]'') dr 



oo 

t 



Im 



53, e 



-iT{A+K) 



(pe-^V)) dr + lm j 

J oo 



4,e 



-iT{A-K) ( - 



dr 



We deal with the first integral, which is equal to (note / G He (A)) 
1 



Im 



-i{A + n) 



-iT(A+K) 



Im ; y (/^3, 
' oo 

1 



1 



-iT{A+K) 



Re 



-Re 



A + K 



i{A + k) 
Pc^n03,^e-*^77 



Re (05, ZT]") - Re / ^05,1 {p/p - ie)z7]'' + z Pc^ \uUF^ - [U, i]er]\ \) dr 

J oo 
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where 



05 = ^^Pc n03. 

We are careful in adding P/^II so that 05 makes sense. We can do so since / G He (A). 
Similarly, the second integral is equal to 

= Re (06, ZT]"") -^^j ('^e, {{P/p- iO)zr]'' + z Pc^ UUF^ - [U, i]^r/ |) dr 



with 



A-K 

We can rewrite their leading terms in the form 



p/n04 . 



Re(05, ZT]'') + Re(06, zr]") = Re j dx {z^5 + #5) iU+ + [/_C)ry 

= Rejdx [{U+ + f/_C) (z05 + z(Pe)] r] 



Re J dx {z(j)s + zcf)^)?]. 
Re (2:07 + z(j)s,T]). 



where we have used (|7.1| ) and (|7^) , with 



f/+05 + t/- 



The remaining integral has the integrand 

= - Re (^05, 1 (p/p - ie)zUr] + z Pe^ C/nF" - [t/, z]^r/ | 
- Re (^06, 1 (^ - ie)zU7] + ^ P/ f/nF« - [f/, i]^r7 | 

= - Re j dx (^yU~"^P05 + /ip06 - ^^(^05 + #6)) 

+ (#5 + #6) ff/nF»- [f/,z]^r7 



- Re y C?X (/i ^J905 + /ip06) t/r] + (-205 + -206) 



where ^277,3 is from F" = i-2"0», + FW, defined in ( pID , 

F« = -lOC - iF{k) - [(ciQ,ImF) + m^]i?n := ^^"0^ + F«"^^^ + F»«^). 
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^,,,3 = -ReJdx {z(j), + z(f)6) Umz^cjyl 

AzrjA = — R-e J dx (^n~^Caz"(l)5 — fJCaZ^cpe) Urj'''^^ 

+ (^05 + #6) [t/nF««3) - C2{z + z)Uir]^^'^] 
^,,5 = -Rejdx {fi-'c^z^<j), - /ic„z"06) f/^^'^ + (/i-'P^^^'^Vs + /ii^^e) ^^r^ 
+ (z05 + #6) [f/nF««^) - {c2{z + z)mr]^^^ + Femr])] 

We first observe that, although e*^ appears in the computation, it does not show up in the 
final results. Also, A^-q^s is a sum of monomials cz^ with = 3. There is no a or r/ in Azri,3- 
Summarizing, we have 

a = {z"^ + z"^) + CzaCiiz + z) + Re{z(j)-j + zcj)^, rj) 

+ I + ^2,3 + A(.a,3) + ^(.„,3)j + [^^'^ + ^2,4 + ^(.a,4) + ^(.„,4)] 

J oo 

We now write 

Apz/^ = + ^2,3 + A^,a,3) + ^(.r?,3) 

= (ciQ, ImAICPC) -4a2oRezc„^°-c,„(z + z) (ciQ,ImF(2)) 
-Re J dx {zcj)^ + #6) f/m^"0i 
From integration by parts, we have 

Apz^ds = apz^- / apz'^ fp{z)ds, ap ^ 



ft 



apz'^ + / ^3,4 + ^3,56/5 
J 00 



where fp is defined in (|6.9| ), and, (/9 = {Po, Pi)) 

^3,4 = -a^2;^(/5o + /5iC)2;-^c„2;" 
^3,5 = -a/3;2^(/5o + /3iC)^-ip(=^'^) 

Note [P] 7^ 0. Also note ap = ap. 

Let a^^-* + a^^-* denote the total of the remaining integrals. We have 



a{t) = a2o{z'^ + z^) + apz^ + c^„a(2; + z) + Re(z07 + #3, ^) + a^^^ + a^^^ (7.3) 
a(^) = / A(^) + ^2,4 + A(,„,4) + ^(.„,4) + ^3,4 ds , (7.4) 

(za,5) 

J 00 



' 00 
ft 

,(5) 
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Now we consider the integrand of a^^\ They are terms of order z^. 
yl2,4 = — 4a2o Re + dibz'^~\ 

A,aA = -c,a{z + z) {ciQ,lm{2\QRaC + 2\Q(r]^^^ + X\C\\)) - 2c,,aRec«z" 
Azr],^ = — Re J dx [n'-^CaZ^cpFj — fiCaz'^cpe) Ur]^'^^ 

+ (^05 + #6) [f/nF««=^) - C2{z + z)f/zr/(2)] 

We substitute a = 020 (-2^ + z"^) + b and t]*^^-* = z^^i^a in the above integrands. Note that 
although terms of order z'^ have the following forms 

4 2 2 2 2 C^) 

z , a , 1] , z a, z 1], arj, zrj^ ' . 

some of them do not occur in the above integrands. We have the following Lemma. 

Lemma 7.1 After the substitution a = 020 (-2^ + z'^) + b and 77^^-' = z°'r]a, the integrand of 
a^^\ can be summed into the form 



B22\z\^ + Re {Aioz^ + A^^z^z + ^2^^} . 

There are no terms of the form 6^, b\z\'^ , or zr]^^\ Moreover, we have 

B22 = |r + O(A^), ^31 = OiX'), A,2 = 0(A) . (7.5) 

Proof. The first part is obtained by direct inspection. Note we deal with 4 in the same 
way we deal with A^^'^^ . The orders of the coefficients are also obtained by direct check, with 
the following table in mind: 

Cq = A, d3o,di2,dos = A^, ^21 = 1, di = 1 
ci = A, R = \-\ 

a20 = A^, a/3 = A^, C^a = A, 07, 08 = A^ 

( z + Xz, aR ~ p^, ~ Ap^ 

We note that most C\z\'^ terms with C = O(A^) are killed by the Im operator. The only 
term that survives, due to resonance, and becomes the main term in i?22; is from the last 
term of A^^^: 



(ciQ, Im AC'r/(2)) = (^c,Q, Im Xz\lz^7]2o) + O(A^) \z\^ + ^ 0(A 



1^1=4,7^(22) 
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where the first term is equal to yFlzl''. We also note that the most dangerous term, = 
z^zu\_ + O(Ap^), in all integrands, only contributes O(A^) to -822- Q.E.D. 

7.3 Integration of 0{z^) terms 

Next we proceed to integrate out those oscillatory terms in a^^\ 

A40Z + A3iZ-'z+Ab2bz ds = —— H \ -— 

—AtK —2lK —ZIK 

where ^(^P^) = bp"^ + 2bpp = 0{z^). Let 

a,o = ^ = 0{\'), a3i = ^ = 0(A3), a,^ = ^ = 0(A) . 

—4tK —ZIK —2tK 

Then 

a^^^ = Re {a^Qz'^ + a^iz^z + ab2bz'^^ + / 5 ds 

J 00 

v44,5 = - Re |a4o2^/4o(^) + 031^:^^/31(2;) + a^^^ipp'^] 

Since Vjf' = f/~^e~*^ e~*"^*f/?7o is larger than 772-5 when time t is of order 1, we also need 
to integrate out terms of order z'^ri^\ as we did for (|6.12|) in the equation for p. Specifically, 
we have terms of the form 



\a\=2 

(for different 0), where 0^2^^ and ^^2^^, 5 are similar to p^ and 5^4, respectively, 
a,^ ~ Re (4., .°e-» ^f^o^ 

\a\=2 ^ J 



|a|=2 ^ J 



(for different </>). Here ^_Q^^j^]^ indicates that we have some resonance effect, (as in p^ and 
(74), and hence these terms has slower dacay in t: 

\a,2^o\ < ii^oii {ty^ {ty^'" , i^2r,o,5i < ii^oii {ty'^^^ {ty^'" ■ 

Now we define 

B5 = B^iz, b, r/(2)^ ^(3)) = ^(5) ^ ^^^^ ^ ^^^^^^^ _^ ^^^^^^^ ^ ^^^^ _ ^^^^^ ^ ^^^^^^^ _ 
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i?5 is of order 0{z^). It is a complicated polynomial in its arguments, which including U, 
U^^ and {{A — ■), but it does not contain e*^. 

We conclude 

a = (^^ + z"^) +b 

b = Czadiz + z) + Re{z(t)^ + ^08, 77) + Re {0402;'^ + 0312;^^ + 052^^} + az^-r^^ 
+ / B22\z\^ + B^ds . 

J 00 

We have obtained the main estimate Theorem assuming the estimate ( [5. 21 ) on h. We 
now need to prove the existence of the solution and check the assumption on h. Define the 
mapping S{h){t) = S'r(6)(t): 



S{b){t) = Czaa{z + z) + Re(2;07 + zcf)^, r]) + Re {0402;'^ + asiz^z + ab2bz'^} + a^^r 

+ [ 522^1^ + 55C?S . (7.6) 

Jt 

Recall the class Bt of b 

BT = {bit):m\<D{t}-\0<t<T} , (7.7) 

where D = 2B22/T = 0(A). Our goal is to show that S{b) maps Bt into itself. More 
precisely, we have the following Lemma. 

Lemma 7.2 Suppose that M{t) < 2 and \b{t)\ < D , < t < T. Recall D = 2B22/T = 
0(A). Then we have 

\s{b){t)\ < c,{D) {t}-''' + ^ {ty' < ^ {t}-' . 

Hence S map Bt to itself. There is a similar statement for a{t). 

Assuming the previous bound on b, we can also estimate 6. Since 6(t) is given by (cf. 



K 

we have 



eU) = '^lmz+ [ -—lm{fi~^p) + Feds 
Jo 1^ 



|e(t)|<01og{t} . 

This estimate will not be used in the rigorous proof, but it provides an idea about its size. 



56 



8 Contraction map 



We review what we have so far. Up to now, we have not shown the existence of the solution 
to the equations setup in section 4 with the boundary condition ( [4. 231) . Notice we cannot 
conclude this from the existence to the Schrodinger equation as the boundary condition of a 
is set at the time t = oo ( [4.23|) . We have obtained the main estimate Theorem |5.1j assuming 
the estimate ( [7.7|) on b. We also proved a bound on the map S in Lemma |7l^ , again, assuming 



an estimate on b. We now need to prove the existence of the solution and prove this bound 
on b. We shall achieve both goals simultaneously by proving the map S" is a contraction 
map on the space B = Bt=oo- Once we have proved this, we have constructed rigorously the 
solution needed in part (1) of Theorem ^Tl] and established all the upper bounds in part (1). 
To conclude part (1) of Theorem |1.1| , it remains to prove the lower bound in (. The size of ( 
is given by z and thus by p. From part (s) of Lemma B]^, we conclude the lower bound of C 



provided that we can check the bound on g. Since g is given explicitly in (6.12) and we can 
estimate all terms in g from the upper bounds in part (1) of Theorem This concludes 



part (1) of Theorem |l.lj assuming that S" is a contraction on the space B. The rest of this 
section is a proof that is a contraction. 

We first recall the setting. For each b & B, we can solve our system to get 

z,p,q = Si{b), T] = S2{b), 9 = Ssib) . 

More specifically, for fixed E and /iq -L Q, our system is 

q = ^21 + ic^bq + g{b, r], 9; t) , 

q = p + 0(2;^), given by ( |6.14| ), z = > 

g(0) given by p{0) = z{0) = (t>. Re ho) + i{u, Im ho) , 

fj{t) = e-'^% + y"*e-*^(*-^) |e''t/nF« - e''[U,t]97]'j ds , 

r^ = U-'e-''fi, r]{0)= P/ /lo , 

9{t) = —lmz+[ Fe{b,r],9-s)ds, Fe = -—lm{fi~'p) + Fe . 

Recall F^ is given by ([4.8|) , and that q = p + cfiz" + cfibz + cfiz^ — p^, for some constants c. 
In this system all a are already replaced by 020 (-2^ + z^) + b. Hence we do not have a. We 
will also write 

p=\q\, q = pe''^. 
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From the equation of q we can write u explicitly as 



uj{t)= f [Im(52ip^ + C46 + Im(^/g)] rfs (8.1) 
Jo 



Assuming the estimate on h in Lemma 7.2, we can bound uo by 



\u{t)\<C\og{t} . 

This provides an idea on the size of tu, but we shall not need it in the following rigorous 
proof for h. 

Return to the proof that S{h) is a contraction. Given 6, h' G i3, we can define t], 6 and 
z', T]', 6' correspondingly. We use 5 to denote the difference of quantities in these two sets. 
For example, 

Sb = b — b' , St] = 1] — T]', 6\z\ = \z\ — \z'\ . 

If 

mt)\<6o{tr' , 

we want to show 

\5S{b){t)\<l5,{ty' . (8.2) 

Let us define 

iV(T)= sup \{ty/'^^\5\z\{t)\ + {tY/' {log {t})-'\5z{t)\ 

0<t<T >- 



.(3) 



+ 



(log {t}r'mt)\} 



loc 



The key here is that 6\z\ has a faster decay than \z\. We first want to show that N{T) < C6q 
uniformly for all T. Our strategy is to show that 

N{T) < Ce^NiT) + C5o (8.3) 

for all T. After this is proved, by choosing Eq sufficiently small we have N{T) < C8q. We 
also have 



\SS{b){t)\ < {t}-'N{t}-''^\og{t} + {t}-"^5,{t}-' 

+ I {sY^'^ N{s}-^'^-'' + {s}-^~''{N + 5^) ds 
J 00 

<C6o {t}-'~^<l5o {t}-' 
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which shows (^.21). Now we focus on proving 



Let t <T. We will write = N{T). Also, when we say things like \S{p'^)\ < C\pSp\, 
what we really means is < C\pSp\ + C\p'5p\. 

1. We first estimate 5g. By (|6.13|) and the definitions of P^'^\ gi, g2, and g^, we have 



\Sg{t)\<C{N + 5o) {t} 



-3/2-cr 



Note C{N + Sq) < Ci if is sufficiently small. Hence 

{ty'^^ < p{t),p'{t) <m{ty'^^ . 



By Lemma |6^ 



\Spit)\<C{N + So)e''{ty^''-''^/' 

2. From the equation for uj we can bound the variation of uj by 

\Sujit) I < r C\5{p') I + C\Sb\ + C\Sg\ {s}'/' + C\g/q'\ \5q\ ds 
Jo 

< [ C{sy'-''e''N + C{sy'6o + C{N + 6o){s} 
Jo 

+ C{N + 5o) {sy'^^-''^' N {s}-^/' log {s} ds 

< [Ce''N + C5o]\og{t} 

Hence the variation of z is bounded by 

\6z{t)\ < \p6u\ + \6p\ < {Ce^'N + C5o){ty'^^ log {t} . 

3. Since 9{t) = C2[{z + z)Yo + /J Fq ds, we have 



3/2-0-+1/2 



\Se{t) \ < C\5z\ + / SFeds < {Ce^'N + C5o) log {t} . 
Jo 

4. To estimate Srj, 



we note 



5 (e'^UUF^ - e'^[U, t]eri] (s) < C\59\ {sy^ + (Ce^N + CSq) {s}' 



hence 



* C 



\m)\\L^< I j^:^ACe''N + CSo){sy'log{s} ds<Ce'^{N + 6o){sy'/'-''' 
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5. To estimate 5^72-5' 



Jo 



Hence 



S4%it) . <{C6''N + C6o){s} 



-1-3(t/4 



loc 



Note that we estimated ||77||^2 in section 5 in order to estimate || |'7p'7||2,8/7, which is used in 



estimating 



~(3) 
^5 



loc 



. However, we do not need ||5?7||2,2 here since 



Wv)\ 



L8/7 



Compare with Subsection 5.3. Other estimates are similar to those in that subsection. 
Since 4% = U-^e-''fjSl,{t), 

Sr^!i\{t) <\5e\ 4%{t) + <Ce'^/\N + So){sr'-/' 



loc 



loc 



5. Summarizing, we have shown (|8.3| ). Choosing Eq sufficiently small we have < CSq. 

6. The above shows that the map S{b) is a contraction mapping, hence there is a fixed 
point b = S{b), which gives us a solution. 

9 Dynamical renormalization of mass 

As explained in the introduction, we expect that the solution to be of the form 

with a changing -E(t), and we use [Qe + clRe + /i]e'®'-*'' as an approximation. In Theorem 
5.1 we conclude that M{T) < 2 assuming \b{t)\ < D{ty' for t G [0,T]. We cannot 
extend this result beyond T. The reason is that, since E{t) is changing and we expect that 
E{t) will converge to some E^o which is likely to be different from E{0), after certain time 
[Qe + clRe + /ije*®*-*-* is no longer a good approximation and \a{t)\ will no longer decay. To 
overcome this difficulty, for each time step T = kAT, we propose to re-choose E = Ek so 
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that the new a{t) satisfies a{T) = 0. Then we prove that the new a{t) has the same estimate 
as the old a{t) if our time increment AT is sufficiently small. Then we estimate the change 
of E{t), especially taking into account of its oscillation, by studying a{t), and prove that 



E{t) does converge. The rest of the proof essentially follows that of Theorem |lT 

9.1 Renormalization lemma 

Lemma 9.1 Suppose Wip — Qe\\2 ^ 1? then we can rewrite 

^ = ((1 + a)QE + k)e'' = {Qe + a Re + h)e'' 
uniquely with h,k Q. Moreover, a, a' , 9, k and h are small. 

Proof. First we can write ip = {1 + a)QE + ki, with ki Qe, and a a complex number. 
Since ip — Qe is small, a and ki are both small. Hence we can find small real a and 6 such 
that (1 + a) = (1 + a)e'^. Let k = kic''^ . Then we have %p = {{I + a)QE + k)e^^ . We 
have Re = c^^cqQe + ^eRe, and hence aQE = ciCQ^a{RE — HeRe)- Let a' = ciCq^q and 
h = k — a' HeRe- We have h ± Qe and gQe + k = a' Re + h. Q.E.D. 



Lemma 9.2 Let tq and 6q be sufficiently small and suppose E E I\ with dist{E,dI\) > 
CtqX. Suppose that ip = [Qe + clRe + h] e*®, with |A^"'^a| = r < tq, \\h\\j^2 = S < 6o. Then 
we can find E with \E — {E + a)\ < rA/2 such that, if we write uniquely 



{9.1] 



using Lemma \9. 1\ , then a = and we have 
{E, 6) is the unique solution of 



h-h 



L2 



< t/2, and |0 - 6| < r/2. In fact. 



^ - Qe' e*®' , Qe 



for 11^ — 11^2 sufficiently small. 



<2Aro), 



(9.2) 



Proof. We will define a sequence \Ek\^^^ 2 s--- which converges to E. Let E\ = E + a. We 
can write uniquely 



= [QE^+aiRE^+hi]e'^\ Hi^Qe^ 



using Lemma |9T|. Denote 6Q = Qe + clRe — Qei and 50 = — ©i. We have \\Qe — QeiW < 
0{aR) < C\a\X-^ < Ct, and ||5Q|| < Ca\dlQ) < Ca^X'^ < Ct\ Note we have 



aiRE, + h, = QE^ie'^"^ - 1) + [SQ + h] e 
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Hence = lm{QEi, ai-R^i + hi) = 0(56) + 0{{Qei, SQ + h)), and hence 

|5e| < 0{6Q) + C\{Qe, h)\ + C\{Qe, - Qe, h)\<CT^ + + Ct6 



Also, 



hence 



Finally, 



O(A-i) ai = {QE,,aiRE, + h) = c^\e'''' - 1) + 0{iQE,,5Q + h)) , 



A">i| < C\e'^® -1\+Ct' + Ct5< Ct' + Ct8 . 



\\hi-h\\^.2 < \\hi- he'^®\\^, + \\he'^® -h\\^, 

= ||(Qi?i(e*'® - 1) + W)e^'® - aiRE,\\^, + \\h\\^, |e^^e - 1| 
< Ct^ + CtS . 

If we choose Ct + C6 < 1/3, then we have 

A~"^|ai|, \SQ\, \\hi — h\\^2 < t/3 . 

Now for /c > 2 we define = -Efc-i + Ofc-i, and define a^, and ©jt correspondingly by 
Lemma |9.1| . We follow the previous estimates to get 

Note that the size of may increase in the process, but since 

C\ak\ + C WhuW^^ < C^l^ + + < C|a,._i| + C Wh^-iW^^ < ■ ■ ■ < 1/3 

our condition for estimates is always satisfied. Hence E + at converges to a limit E with 
|E - E| < 3Ar/2 and \E - E - a\ < Ar/2. Writing ip in the form (|]lD with respect to E, 
we have a = lim^ = and |0 — 0| < Ylk ~ — We conclude 6) is a 



solution of (U). 

To show the uniqueness of ( |9.2|) , we first note that it is locally unique by inspecting the 
equations obtained by taking derivatives of ( |9.2|) with respect to E' and 6'. Now suppose 
we may write 

with hi ± Qi and h2 -L Q2, . Then we have [Qi + hi] = [Q2 + h2]e^^^. Since both hi and 
/12 are small, taking projection on Qi we get 6Q is small. The local uniqueness implies the 
claim. 

Q.E.D. 
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9.2 Proof of Theorem [TT^] 

First we proceed an induction argument to find the desired renormalization at each time 
step and the corresponding estimates. Denote the space of initial data 



Y = H'nw 

Our initial datum is ipo = Q^in + «in-R_Bin + ^in, with \\ipQ — QEiJly — C^^e. Applying Lemma 
9A\ to ipo, we can find Eq, h^, and 9o such that ipo = [Qeo + ^o] Moreover, we have 



1^0 - ^inl < C-^\ain\ < C-^Xe, and hence ||(5in - (5oe'®°^°^ ||y < ^e. This is the beginning 
of our induction argument. We will choose AT > sufficiently small. Then we define 
Tfc = kAT and will define Ej., hk, etc. at time T^. Also, we abbreviate Qk = Qe^, Rk = Re^^ 
etc.. 

Note, Theorem 5.1 remains valid if we replace the assumption \b{t)\ < D {t}^^ by \a(t)\ < 
D {ty^. Also recall that D = 0(A) by Lemma ffl| . 

Assume for = kAT we have found Ek, ak{t), hk{t), Qkit) so that 

f ^{t) = [Qk + ak{t)Rk + hk] e*®S EkElx, 

hk{t)±Qk fortG[0,Tfc]; \\^o - Qke'''''^^'>\\y < e , (9.3) 

. \ak{t)\ <D{t}-\ afc(Tfc) = . 
This is our induction hypothesis. By Theorem 5.1 we have Mk{Tk) < 2. Since ak{Tk) = 0, 



we have 



ak{t) = [a2o{Ek){zl + + " " "It + / 522(^fc)kfc|^ + 



(9.4) 



The direct estimate of Lemma |7.2| gives us 



\akit)\ < lD{t}-' for t G [0,Tfc] . 



Thus |afc(0)| < ^De"^ and we have 



<^£ + A-Vfc(0)|<^£. 

By continuity, there is a cr > such that 

\ak{t)\ < D {ty' for t e [0, Tfc + d] . 

Theorem 5.1 ensures that 

Mk{Tk + a) < 2 . 
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We now return to ( |9.4| ). Since the derivatives of the terms outside the integral (say, zl + z^) 
are bounded by {t}~^, for t G [T^, + a] we have 

\ak{t)\ < CT-\t - T) + CT-^{t - T) (9.5) 

< c,T-\t -t)<\d {t}-^ < \d {ty' 

Z o 

if {t-T)< D/{2C3) and (t - T) < {T} /4 < e~^/A. Now we define 

AT = min {D/{2C3) , 8-^4} . 
Recall Tk+i = + AT. We have thus proved that a > AT and 

\ak{t)\ < lD{t}-\ for t G [0,Tk+i] 



and Mfc(Tfc_|_i) < 2. By Lemma p.2|, we can find Ek+i with 

< ||«fc(Tfc+i)| , (9.6) 

such that we can rewrite 

V^(t) = [Qk+1 + flfe+i (t)i?fc+i + hk+i (t)] e'®'=+i W , hk+i (t) ± Qfc+i 

with afc+i(Tfc+i) = 0. Now we compare these two sets of data with respect to and E^+i, 
for t G [0, Tfc+i]. As in the proof of lemma we have 



<o + cx {T}-^ {ty^^^ < cx {ty^^^ (9.7) 

and since 

ak+i{t)Rk+i + hk+i{t) = [Qk + ak{t)Rk + hk{t)] e^e'^W-^^'^+iW _ Q^^^ 
Taking inner product with Qk+i we get 

A->,+i(t) - au{t)\ < C\Qk{t) - Qk+i\ + C\ak{t)\' + CX {Ty' ||/i,(t)|| < CX {ty'^' 



Hence 



\ak+i{t)\ < \akit)\ + \ak+i{t) - ak{t)\ < {t} ' + {t} ' = {t} ' 

o o 4 



if £ is sufficiently small. By ( p.6| ) and ( |9.?| ), we have 
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Also note 

1^0 - < I |afc+i(0)| < 2De^ < CXe^ . 

Hence Ek+i E I\. We have thus proved the induction hypothesis ( |9.3|) for t = T^+i. The 
induction argument is completed. 

Next we proceed to show that has a limit. Let n = 2^ and T = T„. Since an{T) = 0, 
we have 

rt 

a„(t) = [a2o{En){zl + z^) + . . .] + 522(K)kn|' + ■ ■ ■ 
for t G [0,T„]. In particular, 

|a„(r/2)| <C{T}-' + / 1^221 {sy^ ds < 2D {Ty' 

Jt/2 

(Note that the argument here differs from ( |9.5| ).) Hence, by considering ip = ip{T /2)^ E = E^ 



and E = En/2 in Lemma 9^, we have 



<2|a„(T/2)| <4Z}{T} 



-1 



This estimate shows that the sequence E2k converges to a limit E^o- Moreover, choose ki so 
that 2T2''^AT > e'^ and we have 

\Eo - E^\ < \Eo -Ek,\ + \Ek, - E^\ < CXs^ + CD J] {2^AT}"' < CXe^ . 

k>ki 

Write Qoc = Qe^, Roo = Re^, and 

^(t) = [Qoo + a^{t)R^ + h^{t)] e'®-W 

and compare this set of data with data at T^, t G [0, T^]. As before we get \aoo{t) — ak{t)\ < 
CD {ty^^'^ and hence 

M)\<^-D{ty' . 

Since k is arbitrary, this estimate is true for all t G [0,oo). Since \Eoo — Eq\ < O^Xe^), we 
have — QooWy — ^- Theorem 5.1 shows that 

Moo(T)<2 for all T. 

The first part of Theorem |T]^ is thus proved. 
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Suppose the assumption of part (2) of Theorem [T]^ holds, that is, 



Let ttoo = 0,00(0), 2:00 = -2oo(0), and r^oo = ?7oo(0). Apply Lemma |9J to {^/j, E) = {^po, E-,^) and 
(■?/', E) = {ipo, Eoo) respectively. In both cases, E = Eq. We have 

la-ml < CXe^, |aoo| < CXe^ 
\Eo — -EinI < C|ain|, \Eq — Eoo\ < C\aoo\ 



in "'00 1 



Therefore Qin — Qoo is of order A \a 
Poo = -Pc(-Eoo)noo and Pin = Pc(-£^in)nin, we havc 



0(5^), and hence so is hi^ — h^o- Denote 



||-foo^oo||y ^ ||-Pin^in||y ~\~ ||-Pin(^co ^in)||y ~l" ||(-Poo -Pin)^oo||y 

By the last statement of Theorem 5.1, we have |-2oo(i)| > c {t} for all time t. The last 

Q.E.D. 



statement of Theorem 1.2 is proved 



10 Radiation dominated solutions 



In this section we prove part (2) of Theorem |1.1| . For a given profile ^00, we want to con- 
struction solutions of the form ( |1.6| ) such that 

x(x,t) = /i(x,t)e^®W-^e^*^eoo 
as t goes to infinity. Let W denote the wave operator for C, 

W -.L^ ^ Hc(/:), W(t> = hm e"*^ e'^^-^-^^V • 

We have 



t—*oo 



i{A+E)t 



X, as t ^ 00, 



for general x ^ L"^- See subsection 3.2. For the profile Xoo given in part (2) of Theorem |L1 
we let 



We recall that we set 



^oo = UWXc 



ilj = {Q + aR + h) e^[-^*+^Wl 
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and h satisfies 

dth = -dR - aiQ + C^^'^h - iF{aR + h) - ie{Q + aR + h) . 
To ensure h{t) ± Q all the time, we set 

d= (ciQ,ImF(A;)) , Ci^{Q,R)-\ 

^ = -(1 + cocr'a)-' [a + (cog,Ag'(/i + /i)) + (coQ,ReF(A;))] . 

(Recall Co = {Q,Q)~^-) Then the equation of h is 

dth = Ch + nFaii , 

Ch = -i{Hh + UXQ^U{h + h)} , 

Fall = -iOh -iF - {ciQ, Im F) + iaO^ Ru , Ru := Hi? , 
F = F(aR + h) 

= \Q{2\h\^ + h^) + 2\aQR{2h + h) + 3Xa'^QR^ + X{aR + h)\aR + h) . 

We first rewrite the equation of h. The equation for h is 

dth ^Ch- iOh + HF" , 
F« = -iF{k) - [{ciQ,lmF)+iae]Ru ■ 

Let 

h'' ^Uh . 

Since £ = U~^{—iA)U , we have 

= -lAh'' - meU-^h'' + UUF^ 

= -iAh" - - [u, ijeu-^h" + uuF^ . 

Let h = e*^/i* and use U~^h^ = h, we get 

dth = -iAh + e'^UUF^ - e^^[t/, i]eh . (10.1) 
For a specificed datum infinity, we define a solution by the equations 

h{t)^m+9{t) 

g{t) = f e~'^(*-^) U'UUF^ - e''[U,t]eh\ ds . 

J oo 
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We also let h = U'^e-^^h. 

Suppose we have a solution, then the main term in is of order t~^, hence g ~ 
a ~ t~\ 9 ~ 9 ~ rV2. 

Note that ^{t) is fixed. We consider small in the following norm 

||Coo||i/2n^2,l < s 

with e sufficiently small. Then 

WmWm < CiS, \\m\\w^,o. < C,e\t\-'/\ (10.2) 
Also, for the leading term A|^p(^ in F, we have, for t > 1, 

\\mL^<c\\m\t\\m\\2<ce'{tr\ 

< C{{t)-^'^f • 1 + C {t)-^'^ ■ {{ty^'y < Ce^ (t)-' . 
It i < 1, we simply have HKPCLa < C < Ce^. We conclude 

\\mL.<Ce'{t)-'. (10.3) 
Similarly, if ||5'(i)||^2 < Ce (t)^, one can prove, for example, 

\\\^ + g\\^ + g)\\^,<Ce^t)-\ 1^9^^ < Ce' {t)-'/\ etc. 

Now we proceed to construct a solution. For convenience, we introduce a new variable 
uj = 9. (It should not be confused with the u used in Section 8.) We will define a Cauchy 
sequence on the space 



<^(u;,9,a,g) : [0, oo) ^ R x R x R x i/^ 

\u;it)\ < e'/' {tr'/\ \9it)\ < e'/' {t)-'/\ |a(i)| < £ (^)-^ |b(i)||^2 < £ (t)" 
Here H'^ = IF^'^(M^). We define our map by 

uj^(t) := -(1 + coc^'a)-' [a + (cqQ, \Q\h + h)) + (cqQ, ReF)] 
9^{t) :^ f cuds 

J CO 

a^{t) := f {ciQ,lmF)ds 

J oo 

g^(t) := f e-^^(*-^) {e'^UU{-iF - [(ciQ, ImF) + ma;]i?n} - e'^[U,i]ujh} ds , 

J oo 
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where 

F = F{aR + h), hit) ■= U-\-'' {e~'^'U + git)) . 
Our initial data are 

uj{t) = 0, e{t) = 0, a{t) = 0, g{t) = . 
Given {u!,9,a,g) G A, using this assumption and (|10.2|) -( p!0.3|) , we have 



\\g\'g\\H^<\\9fH^<Cs'(^y 

-2 /.\-3 



ii^ii^.. < {ty 



oo 

t 



\a^{t)\ < / Ce^{sy^ ds<e{t)-^ 

J oo 
J oo 

We have shown that (cu^, 9^ , , g^) G A, that is, our mapping maps A into itself. 

Next we show that it is a contraction. Given [ui, 6i,ai,gi), {002, ^2? ^2) ^72) ^ A, we denote 

So = snp{{tf \5uit)\' + {t) \59{t)\' + {tf |5a(t)| + {tf \\5git)\\^,} 

t 

we know Sq < Se. Notice that —i\^\''^ is cancelled in 6F. 

\\m'g)\\H^<\\9\\l^ mH^<Ce'5o {ty' 

\\5F\\^,<Ce5o {ty'^'<CE5o {ty' 

\5Lu^it)\<cSo {ty'/'<l5'J' {ty"' 

-A/2 ,,\~l/2 



\5e^{t)\< jjl" {sy'" ds<\5T {ty 

\Sa^{t)\< f CeSo{sy'Us<l5o{ty' 
||5^?^WL. < / Ce5o{sy' + Ce'/'5,{sy' ds<\5o{ty' 
Therefore we have 

sup {(t)=' ISco'^it)]' + {t) |50^(t)p + {tf |5a^(t)| + {ty < . 
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These show that our map is a contraction mapping. We conclude that we do have solutions 
h with the main profile e*^*^oo- 
Now 

t) = e'^'h = e*("^*+^)[/-ie-*^/i = e"^^* {U-' + [U, e''] e~''} (e-^^^^oo + g) 
Since \\g{t)\\L2 = 0{t-^) and [U,e'^] = 0{e{t)) = 0{t-^/^), we have 

t) e-'^'U-^e-''^^^, as t ^ oo. 

However, 



, p-iEt i{Il+E)t _ iAt 



oo 



as t ^ oo. Hence part (2) of Theorem is proved 



Remark on radiation dominated solutions to Klein-Gordon equations 

We now sketch a construction for radiation dominated solutions to Klein-Gordon equa- 
tions. We follow the notation in the introduction. For a specified profile r]±, let 

u = ^ + g 

where 

and g denotes the rest. Then we have 

{d^ + B^)i = 0, 

{d^ + B')g = X{^ + gf . 

Hence g{t) satisfies 

g{t) = I' {e^^(*-) - e-^(*-) } ^A(e + g)'ds , 
Since the main source term is 

\\ei<mLuh<cr' 

we get \\g\\2 < iL ^ '^'^ ~ Ct'"^. Then we proceed as in section 10 to construct one such 
solution by a contraction mapping argument. 
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